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Abstract 


An  analytic  method  is  presented  for  estimating  all  the  added  mass  terms  of  a  deeply  submerged, 
submarine  like,  rigid  body.  This  body  may  consist  of  any  number  of  components  (hull,  sail,  fins, 
etc.).  Each  component  is  represented  by  an  ellipsoid  with  three  independently  sized  principal 
axes;  this  allows  the  added  masses  to  be  calculated  analytically.  Ellipsoid  geometry,  orientation, 
and  relative  location  are  chosen  so  that  both  added  masses  and  added  moments  of  inertia  are 
optimally  modelled.  Interference  effects  between  the  main  bull  component  and  an  appendage 
are  approximately  accounted  for  by  using  the  flow  field  around  a  replacement  ellipsoid  for 
the  hull  to  modify  the  flow  at  ths  appendage;  interference  effects  between  appendages  are 
neglected.  The  analysis  uses  incompressible  potential  flow  theory.  It  does  not  account  for  any 
circulation  in  the  flow. 

Calculations  earned  out  using  this  method  are  very  fast.  They  show  that  both  appendage 
thickness  and  hull  interference  can  appreciably  affect  those  added  mass  coefficients  which 
contribute  to  coupling. 


Resume 


On  prlsente  une  mtthode  analytique  pour  estimer  tous  les  termes  de  masse  ajout£e  dans 
un  corps  rigide,  ressemblant  a  un  sous-marin,  immerge  en  profondeur.  Ce  corps  peut  comporter 
un  nombre  quelconque  de  composantes  (coque,  riosque,  ailerons,  etc.).  Chaque  composante 
est  represents  par  une  ellipsoide  posstdant  trois  axes  principaux  de  longueurs  indtpendantes, 
ce  qui  permet  de  calculer  les  masses  ajouttes  de  fagon  analytique.  La  forme,  l’orientation  et 
la  position  relative  des  ellipsoides  sont  choisies  de  fagon  que  les  masses  et  moments  d’inertie 
ajoutes  sont  modelises  de  fa^on  optimale.  L ’interference  entre  la  composante  principale,  la 
coque,  et  un  ajout  est  appreciate  par  l’application  d’un  champ  d’tcoulement  autour  d’une 
ellipsoide  de  remplacement  de  la  coque  qui  modifie  l’tcoulement  autour  de  l’ajout;  1 ’interference 
entre  des  ajouts  n’est  pas  pris  en  compte.  L’analyse  s’appuie  sur  la  theorie  de  l’tcoulement  de 
potentiel  incompressible  qui  n’admet  aucune  circulation  dans  I’tcoulement. 

Les  calculs  effectuts  par  cette  mtthode  sont  tres  rapides.  Us  rtvtlent  que  l’epaisseur  de 
l’jout  et  1 ’interference  de  la  coque  peuvent  influer  beaucoup  sur  les  coefficients  de  masse  ajoutee 
qui  contribuent  au  couplage. 
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A 

B(x) 
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cl>  c2> • • • >  CS 

E,E' 

E(u) 
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H(x) 


H(x) 

U.k 

k,k' 


half  the  lengths  of  an  ellipsoid’s  principal  axes  oriented  along  the  x\  y',  z' 
axes,  respectively. 

the  matrix  transforming  a  vector  description  given  relative  to  inertial 
*o»Vo>*o  11X68  to  description  in  body  fixed  x,y ,z  axe?. 

the  breadth  of  a  hull  component  as  a  function  of  longitudinal  position; 
this  profile  is  assumed  to  be  symmetrical  about  its  centerline. 

the  matrix  transforming  a  vector  description  given  relative  to  body  fixed 
x,y,z  axes  to  its  description  relative  to  the  ellipsoidal  principal  axes, 

s',/,*'. 

appendage  chord  lengths  measured  parallel  to  the  x'  axis;  see  Figure  2. 
complete  elliptic  integrals  of  the  second  kind, 
incomplete  elliptic  integral  of  the  second  kind. 

interference  matrix  containing  the  interference  velocities  at  a  secondary 
component. 

the  height  of  a  hull  component  as  a  function  of  longitudinal  position;  this 
component  is  not  assumed  to  be  symmetrical  about  its  centerline. 

the  z  coordinate  of  the  midpoint  of  the  local  height. 

unit  vectors  associated  with  the  x,  y ,  z  body  fixed  axes. 

modulus  and  complementary  modulus  of  elliptic  functions. 

rolling  moment  and  dimensionless  rolling  moment  on  vehicle;  moment  is 
about  body  fixed  x  axis. 
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M 

\pU*l* 


added  mass  coefficients  giving  the  linear  variation  of  rolling  moment  with 
u,u,u>. 

added  mass  coefficients  nondimensionlized  by  \pl*\  eg,  K £  =  K^/^pi*. 

added  mass  coefficients  giving  the  linear  variation  of  rolling  moment  with 
p,q,r. 

added  mass  coefficients  nondimensionalized  by  \ptl ;  eg,  K £  =  K^/\plh. 
complete  elliptic  integrals  of  the  first  kind, 
length  of  the  vehicle. 

pitching  moment  and  dimensionless  pitching  moment  on  vehicle;  moment 
is  about  body  fixed  y  axis. 
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added  mass  coefficients  giving  the  linear  variation  of  pitching  moment 
with  u,v,w. 

added  mass  coefficients  nondimensionlised  by  \pl* ;  eg, 

added  mass  coefficients  giving  the  linear  variation  of  pitching  moment 
with  p,f,r. 

added  mass  coefficients  nondimensionlised  by  \ptl ;  eg,  =  M^/\ptl. 

matrix  transforming  velocities  and  angular  velocities  specified  relative 
to  body  axes  to  a  representation  relative  to  ellipsoid  principal  axes,  see 
equation  17. 

yawing  moment  and  dimensionless  yawing  moment  on  vehicle;  moment 
is  about  body  fixed  x  axis. 

added  mass  coefficients  giving  the  linear  variation  of  yawing  moment  with 
u,v>. 

added  mass  coefficients  nondimensionlised  by  \pt*\  eg,  =  N^/^pl*. 
added  mass  coefficients  giving  the  linear  variation  of  yawing  moment  with 

P>i>r- 

added  mass  coefficients  nondimensionlized  by  \pl* ;  eg,  =  N^/\plz . 

vehicle  angular  velocities  resolved  about  the  b>dy  fixed  x,y ,z  axes, 
respectively. 

body  component  angular  velocities  resolved  about  its  replacement  ellip¬ 
soid’s  x'ty',x'  axes,  respectively. 

theta  function  parameters;  purely  functions  of  k  and  k' . 

a  generalized  coordinate  in  the  Lagrangi&n  equations  of  motion. 

the  force  associated  with  the  generalized  coordinate  q '  •  in  the  Lagrangian 
equations  of  motion. 

displacement  vector,  using  a  body  axes  representation,  giving  the  location 
of  the  centroid  of  a  component’s  replacement  ellipsoid:  xi  +  yj  -f  *k. 

thickness  to  chord  ratio  for  an  appendage. 

inertia  coefficients  determining  the  kinetic  energy  of  an  ellipsoid  moving 
through  an  ideal  fluid  with  six  degrees  of  freedom;  is  associated  with 
translation  along  the  x'  axis,  with  translation  along  the  y'  axis,  . . . , 
te9  with  rotation  about  the  z'  axis. 

inertia  coefficients  as  above,  except  that  t'Vl  is  associated  with  translation 
along  the  longest  of  the  ellipsoid  principal  axes,  tj2  with  translation  along 
the  second  longest  principal  axis,  . . . ,  t'w  with  rotation  about  the  shortest 
principai  axis. 


v 


T 


the  kinetic  energy  in  an  ideal  fluid. 
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translational  velocities  of  the  vehicle  resolved  along  the  body  fixed  x,  y,  x 
axes,  respectively. 

translational  velocities  of  the  centroid  of  a  body  component’s  replacement 
ellipsoid  resolved  along  its  x',y',x'  axes,  respectively. 

argument  of  the  Jacobian  elliptic  functions. 

speed  of  the  vehicle:  \/u *  +  v*  +  to*. 

translational  velocity  vectors  using  inertial,  body  fixed,  and  ellipsoid  prin¬ 
cipal  axes  representations,  respectively. 

volume  of  an  ellipsoid:  |xabc. 

set  of  inertial  axes,  fixed  in  space. 

axes  fixed  in  the  body  of  the  vehicle:  the  x  axis  points  forward,  the 
y  axis  points  to  starboard,  and  the  x  axis  points  through  the  keel;  the 
orientation  of  these  axes  relative  to  the  inertial  axes  is  given  by  the  angles 
V> ,  0 ,  and  <f> . 
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X ,  *'  = 


X 
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body  fixed  cc  ordinates  locating  the  centroid  of  a  replacement  ellipsoid. 

axes  aligned  with  the  principal  axes  of  a  body  component’s  replacement 
ellipsoid,  such  that  their  origins  coincide  with  the  centroid  of  the  ellipsoid; 
their  orientation  relative  the  body  fixed  axes  is  defined  by  the  angles  $ 
and  Cl. 

an  intermediate  set  of  axes  between  the  x,y,z  and  x',y',z'  axes  used 
for  describing  the  geometry  of  an  appendage,  as  shown  in  Figure  2;  their 
orientation  relative  to  the  body  fixed  axes  is  defined  by  the  angle  $ . 

axial  force  and  dimensionless  axial  force  on  the  vehicle;  force  is  in  the 
x  direction. 
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added  mass  coefficients  giving  the  linear  variation  of  axial  force  with 

added  mass  coefficients  nondimensionlized  by  \pt* ;  eg,  X^  =  X±/ \pls. 
added  mass  coefficients  giving  the  linear  variation  of  axial  force  with 

added  mass  coefficients  nondimensionlized  by  \pt* ;  eg,  X^  =  X^/^pt*. 

lateral  force  and  dimensionless  lateral  force  on  the  vehicle;  force  is  in  the 
y  direction. 

added  mass  coefficients  giving  the  linear  variation  of  lateral  force  with 
u,v,tu. 

added  mass  coefficients  nondimensionlized  bv  \pts ;  eg,  YV  =  Y^j \pt*  ■ 
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added  mass  coefficients  giving  the  linear  variation  of  lateral  force  with 

P.9*^ 

added  maaa  coefficients  nonditnenaionliied  by  \pt* ;  eg,  =  Y^/^pt*. 

normal  force  and  dimensionleaa  normal  force  on  the  vehicle;  force  is  in 
the  *  direction. 


,  Z±  added  man  coefficients  giving  the  linear  variation  of  normal  force  with 

u,v,w. 

added  mass  coefficients  nondimensionlised  by  |p£s ;  eg,  =  Z^/Xpl*. 

Zj,Zj,Zf  added  mass  coefficients  giving  the  linear  variation  of  normal  force  with 

P,4,r. 

Zj}Zj,Zf  added  mass  coefficient!  nondimensionlised  by  \pt* ;  eg,  Zp  —  Z^/^pl*. 


a0,/90)7o  constants  used  in  t)  velocity  potentials  for  flow  around  an  ellipsoid. 

X,n,v  confocal  ellipsoidal  coordinates;  A  =  constant  defines  the  surface  of  an 

ellipsoid. 

p,  v  an  intermediate  set  of  axes  between  the  x,  y,  x  and  x',y'}z'  axes  used 

in  calculating  the  replacement  ellipsoid  geometry  for  an  appendage;  see 
Figure  2. 

0,  t>  coordinates  giving  the  centroid  of  a  replacement  ellipsoid  using  the  p,  v 

axes  of  Figure  2. 

(,<r,r  the  largest,  second  largest,  and  smallest,  respectively,  of  o,6,c. 

p  fluid  density. 

$,0  the  roll  and  sweepback  angles,  respectively,  of  the  x',y ',c'  axes  relative 

to  the  body  fixed  z,  y,  z  axes. 

<p  velocity  potential  in  the  inertial  reference  frame:  v0  =  -V<p. 

tPn  <Pi>  •  •  •  >  Pt  velocity  potentials  giving  the  potential  flow  around  an  ellipsoid  in  each 
of  the  six  degrees  of  freedom. 

Xi>Xs' •  •  •  >Xe  those  parts  of  the  above  velocity  potentials  which  are  purely  functions 

of  A. 

the  yaw,  pitch,  and  roll  angles,  respectively,  of  the  body  fixed  z,  y,  z  axes 
relative  to  the  inertial  x0 ,  y0 ,  zQ  axes. 

<j0, <•>,<■/  angular  velocity  vectors  using  inertial,  body  fixed,  and  ellipsoid  prin  cipal 

axes  representations,  respectively. 


1  Introduction 


DREA  ia  currently  developing  the  ability  to  model  the  maneuvering  characteristic*  of  sub¬ 
marines.  This  will  enable  the  evaluation  of  proposed  candidates  for  the  Canadian  Submarine 
Acquisition  Program  and,  in  the  longer  term,  provide  the  basis  for  a  more  detailed  analysis 
of  the  dynamic  performance  of  the  new  boats.  The  initial  analysis  package  will  be  based  on  a 
numerical  integration  of  the  six  nonlinear  ordinary  differential  equations  describing  the  motion 
in  six  degrees  of  freedom  of  a  maneuvering  vehicle.  These  equations  are  formulated  using  either 
theoretically  or  experimentally  determined  hydrodynamic  coefficients. 

This  memorandum  presents  a  method  for  theoretically  estimating  one  class  of  hydrody¬ 
namic  coefficient,  the  acceleration  coefficients — otherwise  known  as  the  added  masses.  It  is 
the  background  report  to  Reference  1  and  gives  the  complete  details  of  the  theory.  For  the 
convenience  of  the  reader,  it  includes  everything  presented  in  Reference  1. 

Added  mass  coefficients  (added  mass  is  taken  to  include  added  moments  of  inertia)  re¬ 
ceive  their  name  because  they  can  be  linearly  combined  with  the  true  submarine  masses  in  the 
equations  of  motion  to  form  one  coefficient.  Mathematically,  added  mass  is  just  the  propor¬ 
tionality  constant  relating  the  kinetic  energy  in  the  fluid  surrounding  the  vehicle  to  the  square 
of  the  vehicle's  speed,  in  the  same  way  that  vehicle  mass  relates  vehicle  kinetic  energy  to  speed 
squared.  During  vehicle  accelerations,  both  the  kinetic  energy  of  the  vehicle  and  the  fluid  are 
changed,  their  ratio  remaining  constant  in  an  ideal  fluid.  Therefore,  an  accelerating  vehicle 
must  overcome  the  effective  inertia  of  both  vehicle  mass  and  'added  mass’.  The  added  masses 
for  a  submarine  can  be  equivalent  in  magnitude  to  the  actual  submarine  masses,  and  so  must 
be  accounted  for  properly. 

It  is  generally  acknowledged  that  submarine  added  mass  coefficients,  insofar  as  they  are 
important  in  the  equations  of  motion,  are  not  strongly  dependent  on  viscous  effects  (such 
as  circulation  and  boundary  layer  growth).  Indeed,  many  of  the  experimental  techniques 
currently  in  use  for  determining  the  coefficients  make  this  assumption  since  they  ignore  the 
dependence  of  added  mass  on  both  orientation  of  the  vehicle  to  the  oncoming  flow,  and  time  in 
an  unsteady  flow  (stationary  model  oscillations  assume  this,  for  example).  The  assumption  is 
also  consistent  with  most  modern  formulations  of  the  equations  of  motion  in  which  each  added 
mass  coefficient  is  assumed  to  be  a  constant.  Theoretically,  these  assumptions  are  false  given 
the  vortex  structures  known  to  exist  and  the  lifting  surface  Kutta  conditions  that  must  be  met 
in  the  flow;  however,  while  these  flow  phenomena  frequently  dominate  the  steady  state  flow 
field  resulting  from  a  sudden  or  violent  maneuver,  they  are  often  not  present  at  the  beginning 
of  the  maneuver  when  vehicle  accelerations  are  largtct.  All  things  considered,  a  potential 
flow  analysis,  neglecting  circulation,  is  an  appropriate  method  for  estimating  the  added  mass 
coefficients.  If  it  is  found  necessary  to  account  for  flow  structures  and  circulation  effects  in  a 
particular  application,  then  the  equations  of  motion  will  need  to  be  reformulated  to  allow  the 
added  mass  coefficients  to  vary  with  the  state  variables. 

There  are  various  potential  flow  methods  available  for  predicting  the  added  masses,  includ¬ 
ing  exact-numerical  solutions  of  the  flow  about  the  complete  vehicle  configuration.  However, 
a  simple  approach  resulting  in  quick  calculation  times  is  most  appropriate  in  view  of  the  ap¬ 
proximations  discussed  above. 

A  submarine’s  hull  is  generally  responsible  for  the  greatest  contribution  to  the  added  mass 
coefficients.  Hull  coefficients  are  often  obtained  by  representing  the  hull  with  an  ellipsoid. 
This  is  convenient  because  of  the  availability  of  exact  analytical  potentials,  derived  by  Lamb2, 


Z2S 


^  *  -Si  --S  i 


describing  the  flows  sod  giving  ths  added  m sites  for  sllipsoids  moving  in  six  degrees  of  freedom. 
Imlay*  hss  interpreted  Lamb's  work  in  terms  of  a  symmetrically  finned  prolate  spheriod  (ie, 
an  axisymmetric  submarine  without  a  sail);  he  summarises  the  added  mass  expressions  and 
discusses  their  use  in  the  general  ideal  expressions  for  the  equations  of  motion. 

Humphreys  and  Watkinson4  give  a  comprehensive  summary  of  the  usual  formulae  used 
for  estimating  the  added  masses  of  a  multi-component  vehicle;  they  include  a  good  dirrusaion 
of  how  important  each  term  is  to  the  equations  of  motion  (sensitivity  studies).  Tb  simplify  the 
analysis,  submarine  configurations  are  broken  down  to  their  constituent  components,  such  as 
hull,  sail,  and  fins.  The  added  masses  of  each  component  are  calculated  separately  and  then 
summed  with  those  of  the  other  components  to  give  the  configuration  added  masses.  Hull  added 
masses  are  obtained  by  replacing  the  hull  with  an  axisymmetric  ellipsoid.  Appendage  added 
masses  ire  calculated  semi-smpirically  using  flat  plate  models.  The  Humphreys  and  Watkinson 
sensitivity  studies  show  thsee  procedures  give  adequate  accuracy.  However,  Humphreys  and 
Watkimon  do  not.  deal  at  length  with  asymmetrical  hull  shapes  or  configurations  with  sails,  and 
the  sail,  in  particular,  can  have  large  effects  on  those  coefficients  which  contribute  to  coupling 
between  motions  in  different  degrees  of  freedom.  Aucher1 2  estimates  added  masses  in  a  similar 
fashion.  Neither  of  these  methods  account  for  any  hull  asymmetry.  Instead,  they  model  the 
hull  with  an  axisymmetric  ellipsoid  whose  geometry  is  determined  by  an  average  hull  diaimter. 

Another  approach6  uses  strip  theory  to  estimate  hull  coefficients.  Here,  two-dimensional, 
potential,  cross-flow  predictions  estimate  the  local  added  mass  per  unit  length  at  a  longitudinal 
location;  these  local  contributions  are  then  integrated  over  the  length  of  the  hull  to  give 
its  added  masses.  Although  able  to  account  for  hull  asymmetries,  strip  theory  leads  to  an 
overprediction  of  the  principal  added  mass  coefficients  (relative  to  an  exact  3-D  potential  flow 
solution)  of  from  5  to  15  percent,  depending  on  the  coefficient  and  hull  geometry.  Empirical 
corrections  can  be  applied,  but  they  are  data  base  dependent. 

The  approach  taken  in  the  present  theory  is  to  combine  and  extend  the  above  procedures 
and,  as  far  as  is  possible,  apply  the  id<ial  equations  of  motion  iu  a  manner  that  estimates  all 
the  added  mass  coefficients  of  the  vehicle  as  accurately  as  possible.  For  maximum  flexibility, 
each  body  component  is  represented  by  an  ellipsoid  with  three  independently  sited  principal 
axes.  The  ellipsoid  is  positioned  and  oriented  so  that  both  added  masses  and  added  moments 
of  inertia  are  modelled  properly.  As  with  previous  theories,  the  total  added  masses  of  the 
vehicle  are  taken  to  be  the  sum  of  the  added  masses  of  each  of  the  component  ellipsoids, 
except  that  interference  effects  between  the  hull  and  appendages  are  accounted  for  first.  These 
effects  can  be  appreciable,  changing  an  appendage’s  added  mass  contribution  by  as  much  as 
50  to  100  percent  relative  to  a  calculation  which  assumed  the  appendage  was  isolated  in  the 
flow.  The  present  theory  approximately  accounts  for  this  by  using  the  flow  field  around  a 
replacement  ellipsoid  for  the  huh  to  modify  the  flow  at  the  appendage.  Interference  effects 
between  appendages  are  neglected. 

Special  attention  is  given  to  choosing  the  geometry  of  the  hull  replacement  ellipsoid.  This 
is  done  by 

1)  calculating  strip  theory  added  mass  predictions  from  specified  hull  coordinates,  and 
then 

2)  choosing  ellipsoid  geometry  in  such  a  way  that  strip  theory  predictions  generated 
from  the  ellipsoid  coordinates  agree  as  well  as  possible  with  the  hull  strip  theory 
predictions. 
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Thus,  the  hull  replacement  ellipeoid  model*  thoee  geometrical  characteristic*  of  the  hull  which 
are  most  important  in  determining  added  man.  The  added  masse*  for  the  hull  are  obtained 
by  applying  Lamb’s  theory  to  thi*  ellipeoid;  they  can  be  thought  of  as  strip  theory  predictions 
automatically  (though  approximately)  accounting  for  3-D  effects  . 

Appendages  have  only  a  small  effect  on  the  principal  added  masses  associated  with  trans¬ 
lation  and  rotation  in  the  transverse  directions,  but  can  appreciably  affect  and  even  dominate 
those  coefficients  associated  with  surge  and  roll.  Appendages  are  also  important  to  the  so-called 
'off-diagonal'  added  mass  coefficients  which,  though  smaller  in  magnitude  than  the  diagonal 
terms,  are  potentially  important  as  coupling  mechanisms.  Representing  an  appendage  with 
an  ellipeoid,  while  less  satisfactory  than  the  hull  representation,  is  considered  an  appropriate 
approximation  given  the  small  contribution  the  appendage  will  make  to  the  vehicle's  dominant 
coefficients.  Doing  so  also  allows  appendage  thickness  to  be  accounted  for,  and  this  is  shown 
to  have  appreciable  effects  on  some  of  the  coefficients. 

Significantly,  ellipsoid  potential  flow  solutions,  and  the  added  mass  expressions  they  gen¬ 
erate,  require  the  evaluation  of  elliptic  functions  if  the  ellipsoids  have  three  independently  sized 
principal  axes.  This  explains  the  a;  roximation  usually  made  to  avoid  these  calculations:  by 
replacing  the  hull  with  an  axisymmetric  ellipsoid,  only  elementary  functions  need  be  evaluated. 
For  the  present  theory,  fast,  accurate,  analytically  oriented  routines  aru  derived  for  evaluating 
the  requisite  elliptic  functions. 

A  computer  program  has  been  written  to  calculate  the  predictions  of  the  theory  developed 
in  this  report;  its  use  is  described  in  Appendix  B.  In  Section  10,  the  theory  is  applied  to  a 
generic  submarine  configuration  and  results  are  presented. 

2  Th®  Equations  of  Motion  of  an  Ideal  Fluid 

Although  some  of  the  relationships  developed  in  this  section  have  been  presented  by  other 
authors  (Lamb*,  Imlay*,  Abkowits7),  it  is  felt  that  a  thorough  and  complete  presentation 
collated  in  one  place  is  desirable.  This  also  allows  the  reader  to  more  easily  follow  the  extensions 
that  are  made  to  the  above  theories. 

By  deriving  the  equations  of  motion  of  the  fluid  around  »he  vehicle,  expressions  for  the 
ideal  fluid  forces  exerted  on  the  vehicle  can  be  obtained.  One  begins  by  establishing  coordinate 
systems  and  a  set  of  variables  describing  the  motion. 

Consider  the  unsteady  motion  of  a  vehicle  in  six  degrees  of  freedom  having  arbitrary 
translational  and  rotational  velocities.  In  order  to  describe  the  motion,  it  is  conventional*  to 
use  a  set  of  body  fixed  axes,  x,  y,  t ,  which  move  with  velocities  u,  v,  w  in  the  x,  y,  *  directions 
relative  to  a  set  of  inertial  axes  x0,  y0,  z0 .  The  angular  velocity  components  of  the  body  about 
the  x,y,x  axes  are  p,q,r.  Take  the  x0,y0  axes  to  be  in  the  horisontal  plane  and  the  x0  axis 
to  point  downwards.  For  a  submarine,  the  x  axis  points  forward,  the  y  axis  to  starboard,  and 
the  x  axis  through  the  keel. 

Tb  describe  the  orientation  of  the  body  axes  relative  to  the  inertial  axes,  angles  of  yaw  ( ip ) , 
pitch  (#).  and  roll  (^)  are  defined  in  the  following  way7.  Since  finite  rotations  of  yaw,  pitch, 
and  rail  in  any  order  do  not  result  in  the  same  relative  body  axes  orientation,  the  order  of 
the  rotations  is  an  integral  part  of  the  definitions  of  the  angles.  Assuming  the  body  axes  are 
initially  aligned  with  the  inertial  axes  (however  much  the  origins  are  displaced),  an  arbitrary 
reorientation  relative  to  the  inertial  axes  is  obtained  by 
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1)  yawing  about  the  t  axis  through  an  angle  4 . 

2)  pitching  about  the  y  axis  through  an  angle  ® ,  and 
S)  rolling  about  the  t  axis  through  an  angle  4- 

All  rotations  are  in  the  pcaitive  sense  ss  defined  by  the  Right  Hand  Rule.  Thus,  the  matrix  A 
transforming  a  vector  description  using  x0,y0,i0  axes  to  one  using  body  axes  is: 


(1  0  0  \  /  cos  ®  0  -  sin  ®  \  /  cos  4 

0  coed  sin  #||  ^  1  0  I  |  -  sin  V" 

0  -  sin  d  cos  4  J  \  sin  ®  0  cos  9  J  \  0 

(cosfcoeV*  cos  9  sin  4 

—  sin  4 cos d  Jr  sind tin ® cos 4  coedcoed  +  sin  d  sin®  sin  d 

sin  d  sin  4  +  cos^cos^sin#  -  sin  d coe  4  +  cos d sin ® sin  4 


cos  9  sin  4 

coe  d  coe  d  +  sin 4 sin  9  sin  d 
-  sin  d  cos  d  +  cos  d  sin  9  sin  4 


sin  V*  0  \ 

coa  d  0  1 
0 

-  sin® 
sin  d  cos  5 
cosdcos® 


Therefore,  if  v  and  v0  are  the  vectors  specifying  the  velocity  of  the  vehicle  using  body  axes 
and  inertial  axes  coordinate  systems,  respectively: 


v  =  Avf 


■6) 


=  A-1 


Since  A  is  an  orthogonal  matrix,  A-1  =  AT  (where  AT  is  the  transpose  of  A). 

Similarly,  if  <j  and  u0  are  the  angular  velocity  vectors  giving  the  vehicle’s  rotation,  then: 


u  =  Aun 


and  w0  =  Atw  . 


However,  because  of  their  definitions,  d »  ®  (  and  d  do  not  form  the  orthogonal  components  for 
either  of  these  vectors.  The  expression  for  u  is  obtained  by  summing  the  contributions  from 
d.  ®,  and  4  after  properly  transforming  their  associated  angular  velocity  vectors  to  a  body 
axes  representation: 


—  A^Af  At 


+  A^A# 


p  =  -  sin  ®  d  +  d 

q  =  sin  4  cos  ®  d  +  coa  4  ®  (6) 

r  =  coe 4 coa 9 4  -  sin d®  • 

Now,  since  our  interest  is  only  in  the  overall  forces  on  the  vehicle,  which  are  equal  and 
opposite  to  those  on  the  fluid  surrounding  it,  Lagrange's  equations  of  motion  in  terms  of  the 
total  energy  in  the  flow  field  can  be  used.  This  description  uses  generalised  coordinates  (one 
for  each  degree  of  freedom),  all  of  which  must  be  independent.  The  generalised  coordinates 


used  in  this  problem  will  be  taken  to  be  so  that  u,v,  tu.p.g.r  are  dependent 

variables  related  to  the  generalized  coordinates  via  equations  3  and  6.  If  q 1  •  is  a  generalized 
coordinate,  then  the  six  equations  of  motion  are6: 


±(dT\  _  3T 

dt{dq3  )  d?y  "  Q> 


where  T  is  the  total  kinetic  energy  of  the  fluid,  t  is  time,  and  the  Qy  are  forces  on  the  fluid 
associated  with  the  q3-.  Note  that  Q^Sq^  always  has  units  of  energy,  so  that  when  qy  is  an 
angle,  Q3-  is  a  moment. 

Since  it  is  most  convenient  to  describe  the  vehicle  motion  using  body  axis  velocities 
and  forces,  it  is  desirable  to  write  the  equations  of  motion  in  terms  of  u,  v,w,p,q,r  and 
X,Y,Z,K,M,N  (the  latter  are,  respectively,  standard  notation8  for  the  forces  and  moments 
on  the  vehicle  and  are  associated  with  the  x,  y,  x  axes).  This  is  most  easily  done  by  taking  the 
inertial  and  body  axes  to  be  precisely  coincident  at  the  instant  in  time  of  interest.  Then  the 
Qj  will  he  the  negative  of  the  vehicle  forces  and  moments. 

The  terms  on  the  LHS  of  equation  7  are  expanded  as  follows: 

dT  _  dT  du  dT  dv  dTd to  dT  dp  8T  dg  dT  dr 

dq^  du  dq3  dv  dq y  dw  dq ,•  dp  dq3  dq  dq3  dr  dq y  ^ 

Equations  3  and  6  are  used  to  evaluate  the  partial  u  -'vatives  with  respect  to  the  q  -  at  the 

instant  the  axes  are  coincident,  when  i0,y0,z0)4,9,i//  u,v,w:p,q,r: 


d{u,v,w,p,q,r)  _  dv 

d(*o>Vo>zo)  ' 

The  first  term  of  equation  7  becomes: 


(  0,-w,  t/> 

=  I  w,  0,  — u 
\-v,  u,  0, 


du 

Wm) 


0)  — :,0 
r,  0,0 
-q,  0,0 


£  (  dr\  _  £  (dT  du 

dt  J  dt  y  du  dqj 


dT  dv  dT  dw  dT  dp  dT  dq  dT  dr 
dv  dqj  +  dw  dq3-  +  dp  dq}-  dq  dq3  +  dr  dq3 


with: 


dT  dn\  _  (  d  dT\  dn  dT  (  d  dn 
dn  dq^  J  \dt  dn )  dfy  +  dn  dq^ 


and,  again  at  the  instant  of  interest,  one  gets: 

d(u,v,w)  ^  d(p,  q,r)  Q  dv  _  du 

d(x0)y0,i  0)  ’  3(i0,y0,i0)  d(<f>,$,ri>) 

£  d(u,v,w)  _  £  d{p,<l,r)  _  0 

dt  d{<j>,6,ii>)  dt  S(x0,y0,i0) 


1,0,0 

0,1,0 

0,0,1 


dt  d(x0ly0,i0) 


dt  d{4>,9,xp) 


0,  r,  —q 
-r,  0,  p 
9,~P,  0 
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0,  0,-g 

0,  0,  p 

o,  —p,  o 


Thus,  the  equations  of  motion  given  in  Lamb  are  rederived: 


d 

dT 

dT 

+ 

dT 

dt 

du 

r  dv 

9dw 

d 

dT 

dT 

l 

dT 

dt 

dv 

p~dw 

T 

T  du 

d 

dT 

dT 

i 

dT 

dt 

div 

~qdu 

+ 

Pdv 

d 

dT 

dT 

dT 

dT 

i 

dT 

dt 

dp 

w  — “ — 
dv 

*+*  V 

dw 

~rdq 

9  dr 

d 

dT 

dT 

+  W 

dT 

dT 

i 

dT 

~dt 

dq 

Udu) 

du 

~  P  dr 

*r 

r  dp 

d 

dT 

dT 

+  u 

dT 

dT 

i 

dT 

dt 

dr 

~Vdu 

dv 

q  dp 

T 

Pdq 

=  - X 


-  -Y 


(11) 


Since  the  orientation  of  the  inertial  axes  is  arbitrary,  these  equations  are  in  fac;  true  at  any 
instant  in  time.  An  alternative  way  to  write  them  is: 


■ 

(  0 

— r 

q 

0 

0 

0\ 

- 

/  dT/du \ 

r 

0 

~P 

0 

0 

0 

cT/dv 

Y 

d 

-q 

P 

0 

0 

0 

0 

dT/dw 

Z 

dt  ^ 

0 

■-W 

V 

0 

— r 

9 

dT/dp 

K 

w 

0 

-u 

r 

0 

-v 

dT/dq 

M 

• 

\-v 

u 

0 

-9 

P 

0) 

• 

\  dT/dr  ) 

\NJ 

Equations  11  are  linear  in  T,  so  that  if  T  is  written  as  a  linear  combination  of  the 
contributions  from  each  body  component,  then  the  right  and  left  hand  sides  of  the  equations  can 
also  be  separated  into  component  contributions.  This  is  equivalent  to  independently  applying 
the  equations  to  separate  body  components  and  then  summing  the  resulting  component  force 
predictions  to  give  the  total  force  on  the  body.  However,  finding  the  component  contributions 
to  T  is  no  easy  matter  since  they  must  account  for  the  presence  of  the  other  components 
(interference  effects). 

Equations  11  will  be  separately  applied  to  each  of  the  many  components  comprising  the 
vehicle  body,  with  interference  effects  to  be  determined  later.  First,  though,  the  position  and 
orientation  of  the  replacement  ellipsoid  for  each  component  needs  to  be  specified  relative  to 
the  vehicle  body  axes.  We  begin  by  considering  orientation. 

Hydrodynamic  considerations  will  result  in  one  of  the  ellipsoid’s  principle  axes  always 
being  perpendicular  to  the  body  x  axis,  control  surface  deflections  aside.  This  eliminates  one 
degree  of  freedom,  allowing  the  ellipsoid  to  achieve  its  desired  orientation  relative  to  the  body 
axes  with  only  two  independent  rotations. 

Let  the  x\y\z'  axes  be  the  ellipsoid  principle  axes.  Assuming  these  axes  are  initially 
aligned  with  the  body  axes,  their  desired  reorientation  is  achieved  by 

1)  rotating  about  the  x'  axis  through  an  angle  9 ,  and  then 

2)  rotating  about  the  z'  axis  through  an  angle  O . 
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As  before,  all  rotations  are  in  the  positive  sense  as  defined  by  the  Right  Hand  Rule.  This  choice 
of  transformation  means  that  the  z'  axis  will  always  be  perpendicular  i,o  the  body  x  axis.  The 
matrix  B  transforming  vectors  using  a  body  axes  representation  to  ones  using  ellipsoidal  axes 
is: 


(cos  fi  tinfi  0\  /l  0  0  \  /  cos  fl  sin  fl  cos  9  sin  fl  sin  9  \ 

—  sin  Q  cos  n  ollO  cos  9  sin  $  I  =  I  -  sin  fl  cos  fl  cos  $  cos  fl  sin  9  I . 

0  0  1 J  \0  -sin#  cos 9  J  \  0  -  sin 9  cos 9  J 

(13) 

Besides  these  rotations,  the  ellipsoid  will  have  its  centroid  displaced  a  distance  r  =  *1  + 
fj  +  f k  relative  to  the  origin  of  the  body  axes  (using  body  axis  coordinates),  so  that  the 
velocity  of  the  ellipsoidal  axes’  origin  is: 


(u  +  qz  -  ry 
v  +  ri  -  pz 
w  +  py  -  q£ 


In  the  x\t/,  z'  coordinate  system  one  gets: 


v'  = 


=  B  v 


(14) 


(15) 


as  the  translational  velocity  of  the  origin.  Since  the  vehicle  and  its  components  are  undergoing 
solid  body  rotation,  all  of  their  angular  velocities  are  equal  (i.e.,  &  =  u)  and  so: 


Equations  14  through  16  can  be  summarized  in  the  following  matrix  equation: 


where: 


/u') 

V1 

V 

w' 

p' 

=  M 

w 

p 

9' 

9 

\r'  ) 

\r  J 

(17) 
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With  the  various  coordinate  systems  and  their  associated  variables  now  established,  all 
that  is  required  are  expressions  for  the  kinetic  energy  of  the  flow.  These  are  derived  in  Section  4, 
where  the  general  expression  for  the  kinetic  energy  in  the  Sow  around  an  ellipsoid  is  shown  to 

be: 

T  =  \pV  ( tnu'*  +  taaw'2  +  hs™*  +  *44P/2  +  Usl'*  +  *«er'2)  (18) 

where  p  is  flv'd  density,  V  is  the  volume  of  the  ellipsoid,  and  the  are  the  inertia  coefficients, 
determined  by  the  ellipsoid’s  geometry.  Using  this  expression,  the  column  matrix  on  the  LHS 
of  equation  12  becomes: 


'  dT/du ' 
dT/dv 
dT/dw 
dT/dp 
dT/dq 
.  dT/dr  . 


<  du' /du  dv' /du 
du' /dv  dv'/dv 

\du'/dr  dv'/dr 


dr' /du 
dr' /dv 

dr1 /dr 


'  dT/du ' ' 
dT/dv' 
dT/dw ' 
dT/dp' 
dT/dq' 

.  dT/dr' . 


=  f)V  mt  t 


where: 


ftn  0  0  0  0  0  \ 

0  t22  G  0  0  0 

0  0  t8S  0  0  0 

-1  “  0  0  0  t44  0  o  • 

0  0  0  0  t55  0 

V  0  0  0  0  0 

(Note  that  Mr  ^  M”1 .)  Thus,  the  equations  of  motion  for  one  ellipsoidal  component  of  the 
multi-component  body  are: 


d 

dt  + 


0  -r 
r  0 


MrTM 


(X 

Y 

ll  Z 

pV  K 


Since  T  is  a  diagonal  matrix,  MrTM  is  symmetric.  Actually,  this  equation  is  true  for  any 
arbitrarily  shaped  body  providing  T  is  allowed  to  be  a  fully  populated  matrix;  however,  even 
for  this  general  case,  Lamb  shows  that  T  is  symmetric  and,  therefore,  that  MT  T  M  is  always 
symmetric  in  an  ideal  flow. 

Multiplying  through  the  differential  In  equation  20  one  gets: 


since  Mr  TM  is  constant  for  a  body  with  fixed,  immovable  components.  Therefore,  it  can  be 
identified  as  the  added  mass  matrix: 


(  x± 

x. 

x^ 

X , 

X 4 

Xr  \ 

Yu 

Y< 

Y* 

y* 

Y, 

Y* 

z± 

z* 

z* 

zp 

z* 

z> 

K, 

X* 

*4 

X, 

Mu 

m6 

M, 

m4 

M, 

N4 

*r) 

=  —pV  Mr  T  M .  (22) 


Since  this  matrix  is  symmetric,  X 6  =  Zu,  Y+  =  Zit  etc.  Future  calculations  may 

replace  a  term  from  below  the  diagonal  of  this  matrix  with  its  above  the  diagonal  equivalent. 
Equations  23  give  the  expressions  for  the  components  of  the  matrix,  obtained  by  evaluat¬ 
ing  MTTM. 


X,= 


X, 


— pV  [cos3  fltn  4-  sin3  fl  taa] 

—pV  [cos  fl  sin  il  cos  $  (*n  _  *aa)] 
-pV  [cos  fl  sin  fl  sin  ^(*11  ~  *«)] 
yX^  -  2Xil 


iX±  -  IX<, 


n 


Y,,  = 


y,= 


Ya  - 


Yr 


zA  = 


*p 


Zr 


*P 


Ka  = 


Kr 


m4 


K 


Nr 


sX-  -  gx* 

—pV  [sin3  fl  cos3  $  tn  +  cos3  fl  cos3  $  ta 2  4-  sin3  $  tss] 

— pV[cos$sin$(sin3ntn  4-  cos3  fit  „  -  tS3)] 

Wo  ~  *Y* 

z*6  -  £Y^ 

-  yxv 

-pV  [sin3  fl  sin3  $  tn  4-  cos3  H  sin3  ^  t23  -+-  cos3  $  t  ss] 
yZ^  -  2Y4 
2X4  -  iZ<, 

2Y4  -  yX* 

-pV  [cos3  0 t44  4-  sin3  n  t„]  4-  yZ4  -  2Y ^ 

—pV  [cos  fl  sin  0  cos  9  (t44  -  tBB)j  +  gz4  -  zY4 
— pV[coa  Osin  Osin  $(t44  -  tBB)j  4-  gZf  -  2Y+ 

—pV  [sin3  fl  cos3  $  t44  4-  cos3  fl  cos3  $  t55  4-  sin3  $  te0]  4-  2 X4  —  ±Z4 
—pV  [cos  $  sin  $  (sin3  fl  t44  +  cos3  fl  tBB  -  t6C)j  4-  IX+  -  xZ+ 

-pV  [sin3  fl  sin3  9  tu  +  cos3  fl  sin3  $  t65  +  cos3  $  t66]  4-  £Y+  -  yX+ 


(23) 


The  ‘off-diagonal’  coefficients  mentioned  in  the  Introduction  refer  to  elements  of  Mr  T  M . 


Equations  23  are  only  as  simple  as  they  are  because  of  the  special  geometric  symmetries 
of  the  ellipsoid  used  to  represent  a  body  component.  It  ij  obvious  that  a  more  fully  populated 
T  matrix,  representing  a  more  general  shape,  would  lead  to  many  additional  terms  in  all 
of  equations  23.  In  view  of  this,  equations  23  should  achieve  their  greatest  accuracy  if  a 
body  component  is  replaced  with  its  best  representative  ellipsoid  and  the  geometry  of  this 
replacement  ellipsoid  is  used  to  calculate  the  inertia  coefficients.  This  allows  errors  resulting 
from  ‘fitting’  the  component  to  be  directly  monitored  and.  in  fact,  to  be  used  in  determining 
the  best  fit.  On  the  other  hand,  calculating  the  tit  for  a  single  component  from  different 
ellipsoids,  or  from  non-ellipsoidal  geometries,  in  general  invalidates  equation  18  and  will  lead 
to  errors  of  unknown  magnitude.  These  errors,  though  small  relative  to  the  diagonal  added 
mass  coefficients,  might  appreciably  affect  the  smaller  off-diagonal  terms  which  contribute  to 
coupling. 

Care  in  fitting  the  large  hull  component  of  a  submarine  with  its  replacement  ellipsoid  is 
particularly  important  as,  in  many  cases,  the  second  order  effects  alluded  to  in  the  previous 
paragraph  are  roughly  the  same  size  as  the  effects  of  the  appendages.  Fitting  replacement 
ellipsoids  to  appendages  is  obviously  less  important. 

liquation  21  is  the  unsteady  part  of  the  equations  of  motion.  Ueing  the  notation  defined  in 
equation  22,  and  upon  breaking  equation  20  into  its  component  parts,  one  obtains  equations  24 
(on  the  next  page),  which  are  the  complete  set  of  expressions8  for  the  fluid  forces  and  moments 
on  a  vehicle  (or  one  of  its  components)  moving  through  an  ideal  fluid. 

In  addition  to  the  added  mass  contributions  to  the  forces  (Xuii,  Xcv,  ...,  N+r),  which  are 
only  important  during  unsteady  maneuvers,  equations  24  also  give  the  ideal  forces  associated 
with  steady  state  motions  (the  reader  may  recognize  the  M  equation  term  (X^  -  Z^)uw,  for 
example,  as  a  ‘Munk  moment’10).  However,  many  of  these  ideal  steady  state  terms  are  known 
to  be  strongly  subject  to  viscous  effects  and,  where  possible,  should  be  replaced  with  estimates 
based  on  experiment. 

For  the  complete  body,  equations  24  can  be  simplified  if  the  vehicle  is  symmetric  about 
a  longitudinal  vertical  plane,  as  most  submarines  are.  Then,  9  of  the  15  ‘above  the  diagonal’ 
summed  added  mass  terms  are  zero,  namely: 

X<, ,  Xp ,  X* .  Y* > Y*  >  Zp .  £r ,  K4 ,  M+  . 

Of  course,  their  ‘below  the  diagonal’  counterparts  are  also  zero  and,  in  fact,  it  is  by  considering 
some  of  these  latter  terms  that  one  can  most  easily  see  that  a  term  must  be  zero.  Thus,  for 
conventional  submarine  shapes,  equationr  24  reduce  to  equations  25. 
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X  as  X6u  4-  X*(v  -  ur)  4-  X^ftu  -!•  ug)  4-  X^p  +  X4q  4-  X+r 

-  Y^vr  4-  Y^(vq  -  u it)  -  Yprp  4-  {Zt  -  Y4)qr  -  Y+r*  +  Z+wq  +  Z4 gs  +  Z^pq 

Y  -  X6  (it  +  t>r)  4-  Y^v  4-  Y+  (tu  -  vp )  +  Ypp  +  Y4q  +  Y^r 

4-  X6ur  4-  X*(tur  -  up)  +  (X,*  -  Z+)rp  4-  X4qr  4-  X+r*  -  Z+wp  -  Z^p1  -  Z4pq 

Z  —  X^ (it  -  tug)  4-  Y^(v  4-  wp)  +  Z+w  4  Zrf  4-  Z4q  4-  Ztr 

-  X*ug  +  X6(up  -  vq )  4-  (Y4  -  Xjpq  -  X4q 2  -  X+qr  +  YiVp+  Y*p3  +  Y+rp 

K  =  XjU  4-  Y^(v  -  wp)  4-  2p(«>  4-  vp)  +  Kpp  4-  K4{q  -  rp)  4-  Kf(r  +  pq) 

-  X4UW  +  Z'Auv  -  X4ur  +  X+uq  +  ( Z ^  -  Y^vtu  +  Y-^  (v2  -  tu2) 

-  ( Y4  -  Zt){yr  4-  tug)  +  {Y+  4-  Z4)(vq  -  tur)  -  (M4  -  N+)qr  +  M^(g2  -  r2) 

M  =  X4(u  4-  tug)  4-  Y4v  +  Z4(w  -  uq)  +  K4{p  +  gr)  4-  M4q  4-  M+{f  -  pq) 

+  (Xa  -  Z*) utu  4-  X^vtu  4-  Xtif(wi  -  tt2)  4-  (X^  -  Z*)(tup4-  ur) 

-  (X+  +  Z^){up  -  1 ur)  -  Y^uv  4-  Y^vr  -  Y+vp  \  (X, *  -  N+)pr  -  K+(p2  -  r2) 

N  =  X+(u  -  vr)  4-  Y+(v  4-  ttr)  4-  Zf  tv  4-  Kf(p  -  gr)  4-  Mf(q  4-  pr)  4-  N+r 

-  (X*  -  yjuu  +  X6(«9  -  v2)  -  X*utu  -  (. X „  -  Y4)(uq  4-  vp) 

4-  {X4  4-  Y^)(«p  -  vg)  4-  Y^utu  -  Z^tug  4-  Z4wp  -  ( K ^  -  M4)pq  +  K4(p 2  -  g2) 


X  =  X^tt  4-  X^  {w  4-  ug)  4-  X4q  -  Y"vvr  -  Y^rp  -  Y^r2  4-  ^tug  4-  Z-g2 

Y  =  Y^v  4-  Y^p  4-  Y^r  4-  Xatir  4-  X^(tur  -  up)  4-  X4qr  -  Z^wp  -  Z4pq 

Z  =  X^{ u  -  tug)  4-  Ztftu  4-  Z^g  -  Xuug  -  X4q2  4-  Y^vp  +  Y^p1  +  Y+rp 

K  =  Yj(v-  wp)  4-  XpP  4-  Kf  j  4-  pg)  4-  X^uv  -  X4ur 

4-  (Z^,  -  y6)vu  4-  (y,  4-  Z^)(vg  -  tur)  -  (M4  -  N+)qr 

M  =  X4( u4  tug)  4-  Z4( tu  -  ug)  4-  M4q  4-  (X*  -  Z^)uw 

4-  X*(tu2  -  u2)  4-  Y^vr  -  Y+vp+  (X,  -  X,)pr  -  X,(p2  -  r2) 

X  =  Y^v  4-  ur)  4-  Kt{p  -  qr)  4-  Ntr  -  (X^  -  Y^uv  -  X^vw 
4-  (X^  4-  Y^)(«p  -  vg)  4-  Z4wp  -  ( K ^  -  M4)pq 
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3  Interference  Effects:  Primary  and  Secondary  Components 

The  equations  of  motion,  as  summarized  in  equation  20,  a^nly  only  to  an  isolated  body. 
Thus,  equations  23  do  not  account  for  any  modification  of  the  flow  about  the  ellipsoid  due 
to  the  presence  of  another  body.  In  this  section,  these  interference  effects  will  be  modelled 
approximately  by  superimposing  a  uniform  velocity  over  a  relatively  small  ellipsoid  which  finds 
itself  in  the  presence  of  a  larger  one;  this  interference  velocity  will  be  the  velocity  field  around 
the  large  ellipeoid  evaluated  at  the  location  of  the  centroid  of  the  smaller  one,  without  the 
smaller  ellipsoid  being  present. 

Formally,  a  body  component  will  be  assigned  one  of  two  possible  orders:  either  primary 
or  secondary.  In  general,  one  would  call  the  hull  a  primary  component  and  a  hull  appendage 
a  secondary  component.  Primary  component  added  masses  are  not  corrected  for  interference 
effects.  Secondary  component  added  masses  are  only  corrected  for  the  presence  of  the  primary 
component;  all  other  components  are  ignored. 

Thus,  equations  20  and  23  are  correct  for  all  primary  body  components.  For  secondary 
components,  a  modification  must  oe  made  in  the  development  of  equation  20,  at  the  level  of 
equation  17.  In  this  latter  equation,  the  interference  velocity  must  be  subtracted  from  u,v,w 
to  give  the  effective  velocity  of  the  secondary  ellipsoid  through  the  fluid.  There  is  no  correction 
to  p,q,  r  since  the  potential  flow  about  the  primary  ellipsoid  can  have  no  rotation  in  the  field; 
however,  as  primary  component  body  angular  velocities,  p,  9,  r  will  certainly  contribute  to  the 
corrections  applied  to  u,v,w. 

Now,  denote  the  interference  velocity  due  to  a  primary  component  moving  with  velocity  u 
as  Gutti  +  +  tuutik,  velocity  v  as  . . . ,  and  angular  velocity  r  as 

urri+  vrrj  +  tDrrk.  Here,  the  bar  over  a  variable  indicates  the  velocity  is  that  at  the  location 
or  the  centroid  of  the  secondary  ellipsoid.  These  variables  are  dimensionless  if  their  subscripts 
ai«.  u,  v,  or  w,  they  have  dimensions  of  length  if  their  subscripts  are  p,  q,  or  r.  The  velocity 
vector  on  the  right  hand  side  of  equation  17  can  now  be  rewritten: 
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(26) 


where  F  is  the  interference  matrix.  This  multiplication  by  F  is  equivalent  to  subtracting  the 
interference  velocities  from  v,  equation  14.  Equation  17  becomes: 


/u') 

v' 

V 

w' 

=  MF 

w 

P' 

p 

9' 

9 

V  r'  J 

\r  J 

Equation  20  will  then  be  applicable  to  secondary  components  if  MrTM  is  replaced  with 
FrMrTMF,  which  is  the  new  added  mass  matrix. 
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If,  for  a  secondary  component,  the  elements  of  the  matrix  -pV  Mr  TM  are  denoted  , 
X$ ,  . . . ,  Nf  (the  valves  of  these  elements  are  given  by  equations  23  and  they  would  be  the 
true  added  masses  if  there  were  no  interference),  then  the  secondary  component  added  masses 
are: 

x.  =  (1  -  0.)  K(1  -  «.)  -  XX  -  XX.] 

-  D.[XS(1  -  a.)  -  rx  -  y^o.)  -  o„[*S(i  -  aj  -  r>„  -  *J«J 

*.  =  (!-  aj t-Xja,  +  xj(i  -  «„)  -  JtS»J 

-  «.  +  n°(i  -  0.)  -  vje.]  -  o«  [-*X  +  *20  -  ®.)  -  «S»„]  (28) 

N+  -  N$  -  0r  [2  X°  -  X°Gr  ~X°0r  -  Jf9  €D„] 

-  a,  [*y,°  -  XX  -  y»or  -  rX]  -  «r  [2  z,°  -  x°*r  -  rX  -  zX] 

(The  complete  set  of  added  masses  are  easily  obtained  by  matrix  multiplication.)  Note  that  the 
corrections  for  interference  effects  are  proportional  to  second  order  products  in  the  interference 
velocities.  This  is  because  added  mass  is  a  measure  of  fluid  kinetic  energy,  which  varies  as  the 
square  of  the  fluid  velocity.  Thus,  these  corrections  are  potentially  very  important. 

As  with  the  elements  of  the  matrix  MT  T  M ,  the  elements  of  the  inteference  matrix  F 
are  purely  a  function  of  vehicle  geometry,  as  will  be  shown  in  following  sections. 

It  is  noteworthy  that  interference  effects  change  some  of  the  relationships  between  the 
added  mass  coefficients  presented  in  equations  23.  Consider  the  simple  example  of  an  axisym- 
metric  hull  ellipsoid  with  a  sail  positioned  directly  over  its  center  of  buoyancy.  For  roll  and 
sway  motions,  write  the  fluid  kinetic  energy  due  to  the  sail  as  T  —  +  t44p'2)  with 

p'  =  p  and  *>'  —  v~  tp~Cvv ,  the  last  term  in  the  latter  equation  being  the  interference  velocity 
at  the  sail  due  to  the  presence  of  the  hull;  f  locates  the  sail’s  centroid.  (For  this  example,  it 
is  convenient  to  take  the  sail’s  y'  axis  to  be  parallel  to  the  body  y  axis,  unlike  the  general 
notation  presented  in  Section  8.)  If  k  =  1  -  0U ,  v'  =  kv  ~  Xp;  Cu  will  be  negative,  so  k  >  1 . 
Then,  the  time  dependent  parts  of  the  second  and  fourth  of  equations  1  result  in: 

n  =  "*n°;  Yp,X*  =  -K*Yf; 

Kt  =  K?=-pVt44  +  f3Y?  (29) 

As  expected,  Y0  is  proportional  to  k*  ,  and  is  unaffected  by  interference  effects  since  the 
axisymmetric  hull  can  rotate  without  generating  an  interference  velocity.  However,  note  that 

^  — fY0  end  that  the  *1^’  contribution  to  is  not  -fY^  as  in  equations  23.  Symmetry 
of  the  added  mass  matrix  requires  =  K$ ,  and  this  prevents  the  usual  inter-relationships 
among  the  coefficients  from  occuring. 

In  the  above  example,  k  will  be  less  than  2  (a  value  corresponding  to  an  infinitely  small 
soil  with  its  centroid  located  right  on  the  hull  surface).  For  an  asymmetric  hull,  with  its  majc: 
transverse  axis  aligned  with  the  sail,  «c  will  be  larger  than  it  would  be  for  an  axisymmetric 
hull. 
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4  The  Kinetic  Energy  in  the  Ideel  Flow  Around  en  Ellipsoid 

Id  Section  2,  expressions  for  the  added  mass  of  an  ellipsoid  in  an  ideal  Bow  are  obtained  in  terms 
of  the  flow's  kinetic  energy.  In  this  section,  expressions  for  the  kinetic  energy  are  presented. 
Their  complete  derivation  is  presented  in  Lamb*,  beginning  in  article  111.  This  derivation  is 
most  conveniently  carried  out  using  confocal  ellipsoidal  coordinates.  Both  Lamb  and  Whittaker 
and  Watson11,  chapter  23,  present  the  necessary  details  for  use  of  this  coordinate  system. 

The  kinetic  energy  T  of  the  flow  can  be  obtained  if  the  velocity  potential  <p  is  known  in 
the  inertial  reference  frame  of  Section  2.  y?(*0,  y0,  *0)  must  satisfy  Laplace's  equation  V  V  =  0 
subject  to  boundary  conditions  at  infinity  and  the  ellipsoid  surface.  The  flow  velocity  at  any 
point  is  then  v0  =  — and  the  kinetic  energy  of  the  flow  will  be  given  by: 


--//*' 


The  integration  is  over  the  surface  of  the  ellipsoid;  n  is  the  direction  normal  to  the  surface 
into  the  fluid. 

In  the  inertial  reference  frame,  the  equation  for  an  ellipsoid  with  principal  axes  of  length 
2o ,  26 ,  and  2c  is: 

<31> 

The  equation  of  any  confocal  quadric  to  this  ellipsoid  is: 

*» 

_5s - ] - *2_  J - 12 - 1  =  0  (32) 

where  0  is  the  parameter  of  the  quadric.  This  last  expression  can  be  solved  for  0  for  any 
point  (x0,  Vo>*o)  •  Since  the  solution  is  that  of  a  cubic  polynomial  in  0,  there  are,  in  gen¬ 
eral,  three  quadrics  passing  through  (x0,y0,*0)  •  The  confocal  ellipsoidal  coordinates,  de¬ 
noted  X,n,v,  are  the  three  roots  of  this  cubic.  If  (,<r,r  are  the  largest,  second  largest,  and 
smallest,  respectively,  of  a,6.c,  then  — £*  <  v  <  -<r*  <  \x  <  -r*  <  A  <  oo.  The  quadric 
surfaces  X,fx,v  =  constant  we  ellipsoids,  hyperboloids  of  one  sheet,  and  hyperboloids  of  two 
sheets  respectively.  These  surfaces  are  always  mutually  orthogonal. 

Consider,  now,  the  potential  flow  around  an  ellipsoid  translating  through  a  stationary  fluid 
with  velocity  u  in  the  direction  of  one  of  its  principal  axes.  At  the  instant  in  time  of  interest, 
the  centroid  of  the  ellipsoid  is  precisely  coincident  with  the  origin  of  the  inertial  x0,  y0,  z0  axes 
such  that  its  surface  is  given  by  equation  31  and  u  is  directed  along  the  the  x0  axis.  Lamb 
shows  that  a  solution  to  Laplace’s  equation  of  the  form: 

I Pi  ~  *oXi(*)u  (33) 

is  capable  of  satisfying  the  boundary  conditions  for  this  instant  in  time.  He  gets: 


m  abc  dX 

Xl[A)-2-a0Jx  (o*  +  A)A 


where: 


at0  =  a6c  / 
Jo 


o  (a3  +  A)A 


,  A  =  \/(^r+  A)(63  +  A)(c3  +  A) 


A  simple  cyclic  permutation  of  variables: 


a  -»  b  -*  c  and  i0  -»  y0  -»  f0 


u  — »  v  — »  u> 


in  equations  S3  through  35  gives  the  expressions  for  p,  and  ^>s ,  fur  translation  in  the  directions 
of  the  other  two  principal  axes: 

(36) 

«>»=*oX,W<».  X,W  =  J (37) 

where: 

^  (S8) 


where: 

A  ■•**/.”  (* r?W  ■’«■**'/.  (S8) 

a,  i,  and  c,  then,  are  the  semi-axis  lengths  of  the  ellipsoid  principal  axes  parallel  to  the 
directions  of  u ,  v ,  and  w  respectively. 

Consider,  also,  an  ellipsoid  rotating  in  a  stationary  ideal  fluid  with  angular  velocity  p 
about  one  of  its  principal  axes.  At  the  instant  in  time  of  interest,  its  surface  is  given  by 
equation  31  and  the  aforementioned  principal  axis  coincides  with  the  x0  axis.  Again,  Lamb 
shows  that: 

V>4  =  Vo*oXtWp  (39) 


satisfies  these  boundary  conditions,  where: 


_  abc(b »  -  c»)» _  /*» 

j  2(6*  -  e»)  +  (6*  +  e*)0?o  -  7o)  A 


dA 

(6*  +  A)(c*  -f  A)A 


A  cyclic  permutation  in  these  last  two  equations  of  a,b,c  and  x0,y0,z0  with  p,q,r  gives 
y?5  and  <p9 ,  the  potentials  for  rotation  about  the  y0  and  *0  axes  with  velocities  q  and  r , 
respectively. 

Thus,  at  the  instant  in  time  of  interest,  the  linearity  of  Laplace’s  equation  allows  the 
potential  flow  around  an  ellipsoid  moving  in  6  degrees  of  freedom  to  be  written  as: 

<P  =  <Pl  +  Vi  +  V=>3  +  +  <Ps  +  •  (41) 

Equation  30  then  gives  the  kinetic  energy  in  the  flow: 


«■--* ttlJ 


•  ■  >  «  *  v  —  •  - 

»=li=l 

Because  of  the  geometrical  symmetry  of  an  ellipsoid,  only  the  i  =  j  terms  of  this  summation 
are  nonzero.  One  gets: 


T  =  \pV  (tu«»  +  tjjv*  +  t,su;s  +  t44p*  +  i58?s  +  twr *) 
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where  V  =  jxabc  ia  the  volume  of  the  ellipsoid  ind  the  t,,  ere  the  inertia  coefficients: 


5  2 (3*  -  c»)  +  (3*  +  e*)(/J0  -  7o) ' 

1  (o>  - 


S  2(e*  -  a*)  +  (c*  +  o,)(70  ~  °o)  ’ 


Z _ U _ 

66  5  2(a‘-4*)  +  (o‘  +  i‘)(a0-^0)- 

Although  equation  42  was  derived  for  one  instant  in  time,  when  the  ellipsoid  is  aligned  with 
a  purely  arbitrary  set  of  inertial  axes,  the  equation  applies  generally  since  at  any  other  time 
there  will  be  another  set  of  inertial  ar.es  for  which  an  identical  derivation  could  be  performed. 
Note  that,  although  the  notation  is  different,  equation  42  is  identical  to  equation  18  (in  each 
case  the  velocities  are  relative  to  the  ellipsoid  principal  axes). 

The  integrals  in  equations  34  through  38  and  equation  40  are,  in  general,  elliptic  integrals. 
Their  genera!  solutions  are  presented  in  Section  6,  along  with  simplifications  that  can  be  made 
in  special  cases. 

5  The  Interference  Velocity  Field  Around  a  Primary  Ellipsoid 

The  elements  of  the  interference  matrix  of  equation  27  can  be  found  by  differentiating  the 
velocity  potentials  of  the  previous  section.  At  the  instant  in  time  of  interest,  the  velocity  at 
any  point  in  the  flow  field  around  a  primary  ellipsoid  is  v  =  v0  =  -Vtp  (cf,  eq.  2),  where 
<p  is  given  by  equation  41  and  the  ellipsoid’s  boundary  by  equation  31.  The  differentiation  is 
aided  by  the  following  expression  which  follows  from  relationships  between  the  ellipsoidal  and 
cartesian  coordinates3: 

Vy(A)  =  —  /  — ^2 — i-L.  .-fe-j  +  _fB — (44) 

XK  '  2  dX  Vo*  +  A1+  33  +  AJ  +  c*  +  A  )  '  [  ) 

Here,  x(^)  any  function  of  A;  (x0,  y0,  *0)  is  any  point  in  the  flow  field;  is  a  scale  factor 
relating  lengths  between  the  ellipsoidal  and  cartesian  coordinate  systems: 


k?  =  1  /  —fL,  4-  *  il 

U 

1  V(a*  +  A)  (*’  +  A)  (eSxA)7 

If  x(A)  is  identified  with  the  x,(A)  defined  in  the  previous  section,  then  equation  44  can  be 
rewritten  as: 

with: 

72  -  1  72  - 

(  *o  .  .  Vo  .  . 
{a’  +  i,+  i‘  +  X1+; 

:.  +  Ak)P-  «6> 

c,)1  (47) 

1  (2  -  a0)(«3  +  A)  ’  *  [S 

1(3*  —  c3)  +  (32  +  c3)(^0 

72,,  72,  and  72t,72e  are  obtained  by  cyclically  permutating  the  variables  in  and  724  ,  re¬ 
spectively. 


Now,  if  (*o>  9o>  *o)  *•  the  location  of  the  centroid  of  a  secondary  ellipeoid  (coincident  with 
(t,  ff,  f )  at  the  instant  in  time  of  interest) ,  and  if  these  coordinates  determine  X  through  a 
solution  of  equation  32,  then  the  interference  velocities  for  the  secondary  ellipsoid  are: 


-V^/u  =  +  OJ  +  ©„k  =  -Xi(*)i  -  *o?Xi(S) 

=  0.1  +  oj  +  »,k  =  -xt(*)J  -  0ovXi(*) 

-Vy>3/u>  =  Q*1  +  +  ©»k  =  -Xs(^)k  -  I0Vxs(*) 

-V*>4/p  =  0,1  +  0,J  +  w,k  =  -IoX4(^)J  “  9oX4(*)k  ~  MovX«(*) 

-V<p,/g  =  0, i  +  «fJ  +  ©f k  =  -*oX»(^)i  ~  *oX»(*)k  -  loloVxjA) 

-Vio,/r  =  0rl  +  CrJ  +  ®,k  =  -J0Xe(*)l  “  *oXe(*)J  “  *o9ovXe(*)  • 


(48) 


The  Vxi(A)  on  the  right  hand  sides  of  thwe  equations  are  easily  evaluated  in  terms  of 
elementary  functions  of  a,4,e,*o,0o,*o,  as  per  equation  46.  However,  the  x»(^)  are  elliptic 
integrals.  These  constants  are  evaluated  in  the  next  section.  Elliptic  function  theory  must  be 
used  if  the  primary  component  is  asymmetric;  this  is  not  a  difficult  problem  if  the  routines 
necessary  for  calculating  the  inertia  coefficients  for  this  geometry  are  already  in  place.  An 
alternative  approach  would  be  to  use  the  2-D  potential  Bow  field  around  an  ellipse  to  estimate 
the  transverse  components  of  the  interference  velocities,  which  are  really  the  only  ones  of 
interest.  The  advantage  of  the  3-D  calculation  is  that  end  effects  are  properly  modelled,  for 
example,  in  the  region  of  a  submarine's  tail.  As  will  be  seen,  if  the  primary  component  is 
axisymmetric,  then  the  general  3-D  calculations  can  be  carried  out  using  only  elementary 
functions. 


6  Evaluating  the  Inertia  Coefficient*  and  Interference  Velocities 

The  inertia  coefficients  (equations  43)  are  functions  of  a,  b,  and  c  through  the  elliptic  integrals 
c*o ,  Pq,  and  70  (equations  35  and  38).  These  same  integrals,  but  with  different  limits,  need 
to  be  evaluated  when  calculating  the  Xifi)  which  are  present  in  the  expressions  for  the  inter¬ 
ference  velocities.  In  evaluating  the  integrals,  extensive  use  is  made  of  the  theory  of  elliptic 
functions.  A  brief  introduction  to  these  functions  is  given  in  Appendix  A. 

There  are  basically  three  integrals  to  be  solved: 


h 


dX 

(e*  +  A) A  * 


dX 

(<r>  +  A)A  ’ 


dX 

(r*  +  A)A 


(49) 


with: 

A=  v^  +  AN^  +  AKr’  +  A). 

Here,  £,*,r  are  the  largest,  second  largest,  and  smallest,  respectively,  of  the  semi-axes  a,6,c 
of  the  ellipsoid  principal  axes.  Before  obtaining  the  general  solutions  to  these  integrals,  some 
special  cases  are  examined. 
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A  sphere:  f  =  #  =  r 
In  this  cut: 


in  mu  rniii 

7  ./  =/  =f»  /”  «*A  - »  <* 

1  *  *  *  A  (**  +  A)*'*  3  (C»  +  A)8/*  ‘ 

Thua,  Oq  =  /50  =  70  =  2/3  ao  that  tn  =  *„  =  l„  =  1/2  and  t44  =  l,t  =  <fle  =  0. 

A  prolate  spheroid:  r=r 
In  thia  caae: 


rI  =  W~ - . 

A  (o*  +  A)(e*  +  A)' 


_  -<r»  /  2^  ,  g  n/FTa  -  N/fnr^  \ 

f*-**  Vv/?*Ta  v^Tf  v/^Ta  +  y/o  -  o>  /  ‘ 

Using  the  eccentricity  e  =  \Z(*  -  ff*/(,  one  can  duplicate  the  expression  $iven  by  Lamb1 
for  A  =  0 : 

Also: 

T  _T  _c  1  r _ " _ **'  1  r 

*  3  **  A  (o>  +  A) Vf*  +  A  (r‘+A)v/?+A  2/l- 


An  oblate  spheroid:  js# 


A  (e*  +  A)(r*  +  A)3'1  t*  ~  r* 


e*+A 


2  A  -  e*  .  -x  \ 

•*  V  ) 


as  A  — ►  0, 


where  e  =  >/£*  -  r*/{ .  Also: 


dX _ _ 1 

(£*  +  A)*V^TA  (e*  +  A)v/TO  2 


o  A  ■ 


The  only  other  special  condition  that  needs  to  be  considered  is  when  r  =  0 ,  which  is  of 
practical  importance  if  one  wishes  to  model  flat  plate  control  surfaces  (as  is  sometimes  done). 
Unfortunately,  this  apparent  simplification  complicates  matters  (in  ellipsoidal  coordinates,  at 
least).  Not  only  must  elliptic  functions  still  be  evaluated,  in  general,  but  some  of  the  inertia 
coefficients  of  equations  43  become  singular  —  a  condition  that  can  be  remedied  by  considering 
the  expansion  of  the  for  small  r  and  multiplying  the  result  by  V.  This  process  will  be 
presented  following  the  general  solutions  of  equations  49. 


The  general  cate 

It  is  convenient  to  make  the  following  variable  transformation  in  equations  49,  as  per  Munk13 
(it  is  conventional  to  use  the  vrxiable  u  as  argument  in  the  Jacobian  elliptic  functions;  u  has 
nothing  to  do  with  the  velocity  defined  in  previous  sections) : 

_  r* 

A  =  — 5 - (55) 

sn2  u  v  ' 

where,  as  is  usual,  the  Jacobian  elliptic  function  sn(u,fc)  is  written  as  snu.  k  is  the  modulus 
of  the  elliptic  function.  If: 


_  ?  -  <T3 
”  V  -  r2 


so  that  0  <  k  <  1 , 


r  fU  2  j  T  n  [U  8n3  u  J  r  —  fu  sn2  u 

h  =  C  /  sn*  u  du ,  I2  =  C  — j—  du,  I3  =  C  —r-du 

Jo  Jo  dn  u  J0  cn2  u 


where: 


c=  - 


(e2  -  r2)s/J  ' 

Rearranging  equation  55  and  using  equations  A2,  one  gets: 


sn3  u  = 


£3  —  r3 

e3  +  A  ’ 


2  7 
cn  u  =  - 


r3  +  A 


a2  +  A 


I  ~T  ,  2  v  T 

PTa-  dn  “=  F+T 


It  is  convenient  to  rewrite  the  integrals  of  equations  57  in  terms  of  incomplete  elliptic 
integrals  of  the  first  and  second  kinds,  u  and  E(u)  respectively.  Thus: 

/•*•  j  i 

Jq  sn3ud«=  —  [«-£?(«)]  ~  —  [«0-S(«o)]  as  A  — ♦  0,  u  — ►  u0; 
sn3  u  .  —  snucnu  1 


Here,  Jt'3  =  1  -  jfc2  is  the  complementary  modulus.  The  first  of  equations  50  is  obtained 
by  using  equations  A2  and  A5;  the  last  tv  n  by  using  the  periodic  properties  of  the  elliptic 
functions. 

Numerical  values  for  the  elliptic  integrals  are  calculated  using  the  procedures  given  in 
the  appendix:  k  and  k1  are  known;  equations  All  through  A13  give  the  theta  function 
parameters  q  and  q'  (these  parameters  have  nothing  to  do  with  the  angular  velocity  defined 


in  previous  sections);  equations  A14  and  A15  give  the  complete  elliptic  integrals;  u  and  B(u) 
are  then  obtained  using  one  of  the  following  sets  of  equations,  depending  on  the  value  of  q . 

Iff  <0.09, 

1  dn  u  —  y/k'  .  . 

J~2,dn  «  +  Vki  (  ’ 

and  u  is  obtained  from  equations  A8  and  A19;  equations  A24  and  A25  give: 

S(„)  .  4  +  -  sJL  ,  +  O (1000,“) .  (61) 

v  ’  K  dn  u  +  y/lP  2K  1  +  2 q*  cos  ^4u  v 


If  ?>  0.09, 


/  _  J_  dnu->/fccnu  ,g2> 

2  q1  dn  u  +  \/fc  cn  u 

and  u  is  obtained  from  equation  A23;  equation  A28  gives: 

E(«)  =  „  (l  -  +  """(t’cnu  +  '/ldr.u)  +  g  w  .mh  ^r4u  +  0(I000^ ,->«) . 

\  dnu  +  \/fccnu  2ff'  1  +  2g'4  cosh 

(63) 

Equations  A19  and  61  are  good  as  k  -+  0  since  q  — ►  0  at  the  same  time.  When  k  — ►  1 , 
f  — f  1 ,  which  gives  rise  to  the  requirement  for  the  q'  series  of  equations  A23  and  63.  The 
cross-over  point  occurs  when  the  maximum  error  in  these  last  two  equations  equals  that  in  the 
previous  two.  In  equations  A23  and  63,  the  maximum  error  occurs  as  f'  — ►  l/2g'.  As  can  be 
seen  by  replacing  q  with  0.09  (and,  therefore,  q'  with  0.016)  in  the  error  terms  of  the  above 
equations,  the  maximum  error  for  this  worst  case  situation  will  only  be  O(l0-14).  Although 
this  kind  of  accuracy  is  obviously  not  required  in  the  final  estimates  of  this  report,  it  is  useful 
for  the  calculation  of  equations  43  and  59  very  close  to  the  previously  mentioned  special  case 
extremes,  where  the  equations  become  indeterminate. 

An  elliptic  disk:  r=0,  £  >  a 

Only  A  =  0  will  be  considered  here  since,  realistically,  only  control  surfaces  (ie,  secondary 
components)  will  ever  need  to  be  assigned  zero  thickness,  and  the  flow  field  around  secondary 
components  need  not  be  evaluated. 

For  the  elliptic  disk,  some  of  the  inertia  coefficients  of  equations  43  become  singular.  As 
one  might  expect,  they  are  0(l/r),  so  that  T  (eq.  42)  is  finite  when  the  singular  tit  are 
multiplied  by  the  volume,  V.  The  indeterminancy  in  the  expressions  for  Vt w  can  be  removed 
by  expanding  them  in  powers  of  r .  For  this,  the  small  r  expansions  of  the  I’s,  eqs.  57  and  59 
with  A  =  0,  must  also  be  found.  As  r  — ►  0,  one  gets: 


C~jf+0(rs),  k- 
Then,  making  use  of  the  fact  that  u 


~  (i  _“j  +°(r3)> 

-*  K  and  E(u)  — »  E 


+  0{  r2) 


E  when  A,r 


I 


Thus,  with  reference  to  equations  43,  if: 


*u-2  — //  ‘““2  ‘ 

then,  as  r  — ►  0 ,  the  V reduce  to: 


,  _1  (ga-r  a)2(f8-/2) 

44  “  5  2(<r3-r*)  +  (<r3  +  r3)(/3- J8)’ 


rt'u  =  o,  vt'„  =  o, 

4*fr4(g2  -  **) 


™33  = 


4jr£<ra 


yt>  —  . _ 

44  15[(2fJ  -<r3)E-<r*X]  k 


vt'ss  = 


3E  ’ 

15[(£a-2<7 >)E  +  <t*K]  ’ 


^4  =  o. 


(67) 


When  £  —*  <r,  the  last  of  these  expressions  become  indeterminate  since  then  k  — ►  0  and 
K  —*  E.  This  case  must  be  considered  separately. 

A  circular  disk:  r  =  0 ,  {  =  <r 

Again,  A  =  0  only.  Here  it  is  simplest  to  go  back  to  the  case  of  the  oblate  spheroid  (eqs.  53 
and  54  with  A  =  0)  and  let  r  — ►  0,  whence: 


/,  ~  [2  +  0(r»)]  [1  -  (i  +  0( r*))  (|  +  O(r)) 

”  J»  ~  ii + °<r*)  • 


~2-  j  +  0(rs) 


Thus: 

Vt'u=Vf„  =  0,  Vlis=f£s,  Vli,  =»'«'„=  J|  £',  V<i,  =  0. 
Equations  69  are  also  derived  by  Lamb,  but  using  a  different  method. 


(68) 


(69) 


7  Analytical  Expressions  for  the  Added  Masses  of  a  Long  Slender  Ellipsoid 

A  submarine's  hull  length  tends  to  be  large  relative  to  its  transverse  dimensions.  This  allows 
the  expressions  for  the  inertia  coefficients  of  the  previous  section  to  be  expanded  in  inverse 
powers  of  the  largest  principal  axis  length.  One  begins  by  expanding  equation  56: 


Since  k  — *  1  as  £  — >  oo,  q1  will  be  close  to  zero  and  the  general  case  procedures  for  q  >  0.09 
should  be  followed  (see  equations  62  and  63).  The  next  step  is  to  obtain  the  expansion  for  y/k 
from  equation  70,  followed  by  that  for  2e'Q ,  q' ,  etc.  (eqs.  All  and  A12).  Eventually,  one  can 
show  that,  if  A  =  ln(4o/(c  +  6)) ,  then: 

tjj  •»  “(A  —  l)  +  — j  |a(46cA  +  3c8  —  8cl  +  368)  —  4c8  +  7 cb  —  46aj  +  •  •  • 

[2(c  +  A  -  c " 36] 

+  {*(c  +  b) A  [4c(c  +  b) A  -  7c8  -  12c6  -  968] 

+  17c8  +  25ca6  +  51e68  +  516s  }  +  •  •  • 

*33  =  where  s  t„(o,6,c). 

,  (71) 

*«4  ~  ^  lotc  *  { 1 "  [12(c  + b)3 A  - 15c*  -  26cb  ~  15fc3]  +  —  } 

ti5  ~  —  [6(c  +  6)aA  -  3c8  -  20c6  -  568] 

+  j  |l2(c  +  5)aA  [l26(c  +  6)A  +  c8  —  60c6  —  25b8] 

+  3c4  +  598cs6  +  1463ca68  +  794 cb3  +  177b*  }  +  • .  • 

*ee  =  hs(a>c>b)>  where  t86  =  *55(o,6,c) 

Equations  71  are  particularly  useful  in  showing  how  hull  added  masses  are  analytically 
related  to  hull  geometry.  For  numerical  calculations,  one  can  expect  the  error  in  the  expansions 
to  be  of  the  order  of  the  first  neglected  term.  Although  extensive  calculations  have  not  been 
performed,  the  equations  appear  to  givs  satisfactory  engineering  accuracy  for  l/d  >  6;  good 
accuracy  is  obtained  for  l/d  >  8 . 

Since  the  present  theory  evaluates  the  t w  exactly,  it  does  not  make  direct  use  of  equa¬ 
tions  71. 


8  Replacing  a  Component  with  its  Best  Fit  Ellipsoid 

This  step  is  as  important  as  any  so  far.  Benefits  from  the  analysis  of  the  previous  sections 
will  be  maximised  if  a  submarine  component  is  replaced  by  an  ellipsoid  whose  geometry  and 
position  are  optimally  chosen  (due  consideration  being  given  to  the  approximate  nature  of 
these  estimations).  The  following  analyses  show  the  important  characteristics  this  replacement 
should  model. 

A  hull  component 

A  hull  component  replacement  ellipsoid  will  have  its  largest  principal  axis  (length  2a)  parallel 
to  the  x  axis  of  the  conventional  submarine  body  fixed  axes.  Its  width  and  height  are  deter¬ 
mined  by  the  lengths  of  its  other  two  principal  axes,  2b  and  2c ,  which  are  parallel  to  the  body 
y  and  x  axes,  respectively.  All  hull  components  are  assumed  to  be  symmetric  about  a  vertical 
plane  through  their  centerlines.  This  plane  will  be  parallel  to  the  x-z  plane  formed  by  the 
submarine  body  axes,  but  may  be  displaced  from  it  by  an  amount  y.  Although  y  will  u  '  ally 
be  sero  for  the  main  hull  component,  in  general  it  will  not  be  for,  say,  an  outboard  prS  for 
which  the  'hull  component’  representation  is  appropriate. 

With  fl  =  9  —  0,  equations  23  simplify  to: 


Kp  =  -pVUi  +  ?z*  + 

K4 =  -xyZ^ 

K,  =  -S2Y, 

M^-pVtn  +  g'Xt  +  s'Z* 

M+  =  -y  zXt 

Nf  =  -pVtw  +  x*Yi  + 

Equations  71  show  how  sensitive  these  added  mass  terms  are  to  variations  in  a,  b ,  and  c. 
When  multiplied  by  V  =  Jxafcc,  the  lowest  order  terms  in  the  expressions  for  t33,  ts3,  t55, 
and  t6e  are  independent  of  the  dimeneion  in  ike  direction  of  motion  asst  dated  with  the 
coefficients  Yif  ,  M4 ,  and  N+  respectively;  this  dimension  has  only  a  small  affect  on  these 
added  mass  coefficients.  In  other  words,  modelling  the  characteristics  of  the  maximum  cross- 
sectional  area  perpendicular  to  the  direction  of  motion  is  of  primary  importance  to  these  four 
terms;  modelling  the  volume  is  not  important. 

This  is  consistent  wi'ih  the  expression  for  the  added  mass  per  unit  length  of  an  elliptic 
cylinder  moving  in  a  direction  perpendicular  to  its  longitudinal  axis3: 

*ph\  (73) 

where  2b  is  the  length  of  the  principal  axis  (of  the  cross-sectional  ellipse)  oriented  perpendicular 
to  the  motion.  The  dimension  parallel  to  the  motion  has  absolutely  no  effect  on  this  result. 


X «  —  pYtxl 

x4  =  EX* 

X,  =  -  sx* 

Y*  =  -pVtn 
Yp  =  SY, 

Yf  = 

K  =  -P^ss 
Zp  =  VZw 

Zq  = 
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One  can  use  this  2-D  result  in  a  strip  theory  formulation  of  the  added  masses.  If  B(x)  is  the 
local  breadth  of  the  hull  component,  H(x)  the  local  height  as  seen  in  a  sheer  profile,  H(x)  the 
t  coordinate  (in  conventional  body  fixed  axes)  giving  the  midpoint  of  the  local  height  line, 
and  t  the  length  of  the  component,  then  strip  theory  formulations  of  the  following  added  mass 
terms  make  use  of  the  indicated  integrals: 


1  H*{x)  dx  = 

16  2 

— oc 

3 

J  B7(x)  dx  = 

Yt,Nt: 

j  xH7(x)  dx  : 

■  *16  2 
=  £—ac 

3 

Z<j » •Mui : 

4s  = 

J  xB*(x)  dx  = 

3 

(74) 

Y^K,: 

fW(x)H\x) 

j  ,  16  2 

dx  —  £~ac 

3 

Vhs-  h  =  J(x-i)2Bi(x)dx=—aabi 
VtM:  h  =  J(x  ~  i)2Hi(x)dx  =  jj:asca. 

These  I  integrals  are  easily  calculated  from  hull  coordinates.  The  RHS’s  of  the  equations  give 
the  values  of  the  integrals  for  an  ellipsoid  with  its  centroid  located  at  (£,y,£) .  Thus,  the  best 
fit  hull  replacement  ellipsoid  will  be  one  that  best  satisfies  equations  74.  Note  that  X  and  K 
derivative  characteristics  do  not  contribute  to  this  decision  making  process.  This  is  because 
X  derivatives  are  less  important  to  submarine  maneuvering  characteristics  than  the  cross-flow 
terms,  and  because  Vt44  for  the  hull  component  is  invariably  dominated  by  the  other  terms 
and/or  the  sail’s  contribution  to  K ^ . 

In  the  calculations  that  follow,  y  is  assumed  to  be  a  known  constant,  determined  by  the 
location  of  the  hull  component,  a,  b,  c,  x  and  z  are  all  unknown  and  are  calculated  from 
the  information  provided  about  the  hull  component,  as  summarized  in  the  I  integrals.  The 
problem,  of  course,  is  that  there  are  seven  equations  in  five  unknowns.  This  is  resolved  as 
follows. 

i  can  be  calculated  directly: 

(75) 

x  is  taken  to  be  the  average  of  its  horizontal  and  vertical  plane  predictions: 

*  =  (W*.  +  ««/!  (76) 

and  the  errors  associated  with  this  approximation  are  monitored  with: 

«,  =  (*-  =  -(*  -  (”) 
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Ellipsoid  semi-axis  lengths  a,  b,  and  c  are  determined  by  the  remaining  four  of  equations  74, 
by  minimising  the  sum  of  the  squares  of  the  errors  with  which  the  equations  are  satisfied.  Let: 


<78> 

Then,  if  E  =  e|  +  ej  +  «J  +  ej ,  the  three  equations  determining  a ,  6 ,  and  c  are: 


There  results: 


so  that: 


a  = 


f  J  -  f !  I  CS  — 


as  « 

db  ’ 

o 

1! 

wi « 
(0|‘0 

(79) 

■e »  cs(1  + 

«s)  =  -«s(l  +  «s) 

(80) 

V.Y'4 

Vs/  ’ 

2  37j  1  +  v^-VeA-Vs) 

'  16a  1  +  W(V5)  ' 

(81) 

Equations  75,  76,  and  81  define  the  best  fit  hull  replacement  ellipsoid. 

Now,  consider  the  added  moment  of  inertia  per  unit  length  for  a  rolling  elliptic  cylinder, 
which  is3: 

»p(c3-fc3)J/8.  (82) 

Strip  theory  makes  use  of  this  equation  in  formulating  the  following  integral,  which  is  used  for 
calculating  K^\ 


Vt4i:  I4  =  Jt  [H' »(*)  -  B3(*)ja  dx  =  ^a(c3  -  b3f.  (83) 

Here,  again,  the  RHS  of  this  /ration  is  the  value  of  the  integral  for  the  hull  replacement 
ellipsoid.  The  degree  to  which  equation  83  is  satisfied  is  monitored  by: 

U  ^  ”  b 2)2  "  1 '  fR4) 

What  one  really  wants  to  do  is  monitor  the  errors  in  equations  72  compared  to  the  total 
added  masses  (written  as  (  )T)  for  the  complete  vehicle  configuration.  This  can  be  done  by 
monitoring  the  following  relative  errors: 

=  ejn/(y6)T 

<k  =  + /(*,)T  (85) 

eM  =  (— C5^*5S  +  ~  2%£xIZ^j) / (Mj) ^ 

€ff  =  (— £^Vt9ii  ■+■  +  2 i£sfl/ri)/(iV|1)T  . 
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a)  Antisymmetric  bull. 


b)  MyPim^tric  bull  with  straight  vertical  stem. 

Figure  1  Sheer  plans  for  tW°  standard  hulls  and  sail  with  their  replacement  ellipsoids.  The 
hulls  have  identic*!  breadth  plans. 

Here,  it  is  understood  thut  the  added  masses  in  the  numerators  of  these  expressions  are  those 
f  r  the  hull  component;  also,  except  for  c4 ,  the  f,  are  assumed  to  be  small.  These  expressions 
estimate  the  relative  erf^r*  in  the  overall  added  masses  resulting  from  fitting  an  asymmetric 
hull  component  with  a  singlo  ellipsoid.  For  axisymmetric  hull  components,  the  ti  are  all  zero 
iince  then  equations  74  r^du^e  to  four  equations  in  four  unknowns,  and  these  can  be  satisfied 
exactly. 

Figure  1  shows  the  she**  Plans  °f  two  different  submarine  hulls,  along  with  their  best  fit 
ellipsoids.  The  asymmetry  iq  the  hull  of  Figure  lb  is  confined  to  the  nose  section,  the  front  20 
percent  of  the  hull  length. 

For  the  asymmetric  hull,  e  values  are: 

e3  =  0.0l9>  e,  s=  -0-02U,  e4  =  -0.67,  ex  =  -0.009 
and  for  the  hull  and  sail  combination  shown: 

=  0.016,  ez  =  -0.020, 
eH  ®  ^.0°3 ,  eM  =  0.019  ,  eN  =  -0.019 . 

Note  the  extremely  small  effnc*  th®  very  large  e4  value  has  on  cK ;  this  is  because  of  the  large 
contribution  to  K ^  from  tbn  8  ail-  This  asymmetric  hull  gave  the  largest  e3,es,e4  values  of 
any  hull  tested.  This  ie  most  likely  due  to  the  hull’s  asymmetry  being  all  at  one  end,  while  th«s 
ellipsoid  modelling  it  is  evenly  asymmetric  over  its  entire  length. 

Equations  81  are  responsible  for  these  low  e  values.  An  initial  attempt  at  reducing  equa¬ 
tions  74  to  five  equations  in  five  unknowns  was  made  by  assuming  B(x)/H(x)  —  b/c  for  all  z, 


and  using  /4  to  determine  the  difference  between  b  and  c ;  however,  this  resulted  in  e  values 
as  large  as  0.05,  even  though  e4  was  very  small. 

Section  10  gives  an  example  of  another  hullform. 

An  appendage 

Three  types  of  appendages  will  be  considered,  as  shown  in  Figure  2.  Most  appendages,  and 
certainly  the  sail,  will  need  to  be  represented  by  a  type  I  or  II  appendage.  The  type  III 
appendage  could  be  useful  for  modelling  sailplanes  or,  perhaps,  endplates.  It  should  not  be 
used  to  model  a  pair  of  appendages  separated  by  the  hull  if  interference  effects  are  also  being 
modelled,  since  the  centroid  of  the  replacement  ellipsoid  would  then  lie  within  the  hull  and  the 
interference  velocity  at  this  centroid  would  have  no  meaning. 

The  specifications  required  to  completely  determine  the  geometry  and  location  of  each 
type  of  control  surface  are: 

Types  lie II:  (*x,Vi,*i),  (*a»Va)>  (*a>Va)>  (*4>V4)>  * 

Type  III:  9,  elt  (x4,y4),  c4,  (x5,y5),  c8,  t. 

For  each  type,  »re  the  body  axis  coordinates  of  the  trailing  edge  at  the  root  of  the 

appendage.  This  point  is  also  the  origin  of  the  local  x,y ,z  axes.  The  £  axis  »  always  aligned 
with,  though  usually  displaced  from,  the  body  z  axis.  9  is  the  angle  the  y  axis  must  roll 
away  from  the  body  fixed  y  axis  to  bring  itself  into  the  plane  of  the  appendage.  ,  y* » 0)  > 
«  =  2, 3, 4, 5 ,  are  the  corner  coordinates  of  the  appendages  in  the  local  £,  y,  z  coordinate  system; 
they  determine  Q .  (For  the  type  III  appendage,  root  and  tip  cord  lengths  are  specified  instead 
of  leading  edge  coordinates.)  Together.  9  and  fi  give  the  orientation  of  the  /i,  u  axes,  and  of 
the  z',y',z'  axes.  As  discussed  in  Section  2,  9  and  O  are  zero  when  the  y'  and  x'  axes  are 
aligned  with  the  body  fixed  y  and  x  axes,  respectively;  positive  rotations  are  determined  by 
the  Right  Hand  Rule,  with  a  rotation  in  9  being  applied  before  that  in  fi .  The  thickness  to 
chord  ratio  of  the  appendage  is  designated  by  t . 

The  centroids  of  the  ellipsoids  replacing  these  configurations  are  located  at  (£,?),  as 
discussed  below.  These  points  locate  the  origins  of  the  z',y'  axes,  which  are  aligned  with, 
though  displaced  from,  the  /i,i/  axes. 

Note  that  the  n,  v  axes  defined  in  this  section  have  nothing  to  do  with  the  confocal 
ellipsoidal  coordinates  presented  in  Section  4. 

For  type  I,  ft  is  determined  by  taking  the  v  axis  to  be  parallel  to  the  bisector  of  the  angle 
formed  by  the  leading  and  trailing  edges  of  the  appendage.  The  y'  axis,  also  parallel  to  this 
bisector,  is  allowed  to  be  displaced  from  it  to  partially  account  for  asymmetric  end  geometries. 
This  type  I  representation  can  also  be  used  for  triangular  fins  in  which  the  points  (x3,y3)  and 
(z4,y4)  are  coincident. 

For  type  II,  the  p  axis  is  simply  taken  to  be  the  line  joining  the  leading  and  trailing  edge 
root  coordinates;  that  is,  O  =  tan_1(y3/x3).  Everything  else  is  the  same  as  type  I.  Since  this 
type  II  representation  does  not  require  leading  and  trailing  edges  in  order  to  calculate  fl ,  it 
can  accomodate  all  kinds  of  triangular  fins,  including  those  in  which  the  points  (x3,  y3)  and 
(z3,y3)  or  (z4,y4)  and  (0,0)  are  coincident. 

For  type  III,  the  simplification  of  setting  H  =  0  is  made.  This  type  partially  accounts 
for  asymmetries  in  the  appendages  about  the  n  axis  and  the  line  p  —  c1/2  by  strategically 
locating  the  ellipsoid’s  centroid. 
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Type  1:  fi  determined  by  bisector. 


Type  tl;  0 


Type  HI:  0  =  0. 

Figure  2  Appendage  profiles  in  tiie  plane  l  =%  0 


lb  choose  the  optimum  dimensions  end  position  of  a  replacement  ellipsoid,  one  identifies 
those  characteristics  which  are  most  desirable  to  model.  Unlike  the  hull,  where  simultaneously 
monitoring  motions  in  four  or  all  of  the  degrees  of  freedom  is  appropriate,  an  appendage  pri¬ 
marily  affects  forces  perpendicular  to  its  planform  area.  This  suggests  three  degrees  of  freedom 
nr*  important:  translation  along  the  *'  axis  and  rotations  about  the  z'  and  y'  axes.  However, 
only  the  first  two  of  these  are  of  primary  importance  because  of  the  very  important  role  ap¬ 
pendages  play  in  vehicle  rolling  motions  (this  assumes  ft  is  not  large).  Properly  modelling  the 
thickness  of  the  appendage  is  not  important  when  modelling  these  motions.  Of  course,  having 
determined  ellipeoid  geometry,  the  added  masses  for  all  degrees  of  freedom  will  be  calculated. 

If  e(i>)  is  the  local  chord  length  (perpendicular  to  the  v  axis)  of  a  type  I,  II,  or  III 
appendage  and  ?{v)  is  the  u  coordinate  of  the  midpoint  of  the  local  chord  length,  then  the 
force  and  moment  characteristics  that  a  replacement  ellipsoid  should  model  are: 


Is  6 dV  ~  * 

J mi  s  fs  ?(*/)c*(I/) dv  =  . 


Jyi-  f  vc2{v)dv  =  £>^0*6, 

J  s 

J¥l~  f  i/V(v)d»/=  ffo*6s  +  p*^ol6. 
J  s 


(86) 


These  integrals  model  the  distribution  of  the  square  of  the  chord  length  along  the  span  (5) 
of  the  appendage,  which  equation  73  indicates  to  be  an  important  two-dimensional  added 
mass  characteristic  of  the  cross-flows  resulting  from  the  translation  and  rotation  of  primary 
importance.  The  righthand  sides  of  the  equations  give  the  results  of  the  integrals  for  an  ellipsoid 
With  its  principal  axes  aligned  with  the  z',  y#,  *'  axes  and  ita  centroid  at  (p,  p) . 

The  following  equations  give  the  expressions  that  are  uaed  for  c(u)  and  C(v)  in  the  above 
integrals: 


Types  I  &  II  < 


c(*0  = 


'»> 


a 

l/~Vk 


,  .  c,  v—vk  c. 


c(l/)  =  C*  -  (c*  -  ce)  ,  C{u)  = 


__  Vj-V 


V4-Ve 


2  ’ 

*(„) «  £i  +  JiZ*  *4 . 

v  ’  2  vA  —  i/j  2 


<  V  <  l/j 
Vl  <  V  <  ve 


ve<v<vd 


(87) 


f  cW=cr^, 

Type  III  |  5 

[cM=cl-e-^»t 


e(^)  =  j  + 


2»s  +  c8-Ci 

2ys 

?*i±iiz£iu 

2y< 


ilj  <  ^  <  o 
0  <  v  <  y4 


where,  for  types  I  and  II:  va  and  vh  are  the  minimum  and  maximum  of  ux  and  vt ,  respectively, 
l>£  and  vd  are  the  minimum  and  maximum  of  and  vt ,  respectively,  and  ci  is  the  chord 
length  perpendicular  to  the  v  axis  at  v  =  i/<.  Thus,  the  type  I  and  II  representation  allows 
to  be  greater  than  or  less  than  sero,  and  uA  to  be  greater  than  or  less  than  .  The  type  I 
representation  cannot  be  used  if  it  results  in  vc  being  less  than  vk ;  type  II  must  be  used 
instead. 
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Fbr  a  type  I  or  II  appendage,  the  limits  of  integration  in  equations  86  are  from  to  . 

One  gets: 

Jq  *  +  «#**(*'.“*'»)  +  c»(-  -»'.)]/* 

J*i  =  "  [vK*»«-*i)  +  «*e.^i]  A- 

J»i  =  K(*'<“*'»)(*'*+2*'.+*'*)  +  2c.e*(*'.*-*'*)  +  c*(*',-»'JK+2^+^)]/ja 

+  c*K-O(^+2^+^v.+S^+2i/ky.+»/*)]/S0 . 

(88) 

For  a  type  m  appendage,  the  limits  of  integration  are  from  yt  to  y4 : 

Jo  =  [(ci+cic4+e4)04  -  (c*+«ies+e*)Vs]/S 

JMl  =■  -^^0+  [(ci+2cic4+Sc4)(2*4+c4“ci)V4  ~  (c*  +2e1c5+3e|)(2£5  +  c8 -c^slsj/24 

J»1  =  [(Cl  +  2ClC4  +  8«M  -  («!  +  2clC8  +  5ci)y?]/l2 

J*i  =  [(e*+8cie4+®c4)V4  -  (cJ-h3r1c8-»-6c|)y|3/SO- 

(89) 

Finally,  for  each  type  of  control  surface: 


The  ellipsoid  principal  semi-axis  length  along  the  *'  axis,  c ,  can  be  chosen  without  com¬ 
promising  the  above  more  important  calculations.  If  t  is  the  thickness  to  chord  ratio  of  the 
control  surface,  then  the  following  equation  models  the  distribution  of  the  square  of  the  thick¬ 
ness  perpendicular  to  the  n~v  plane. 

J  [<c(«/)]*  dv  =  ^c*6. 

Here,  again,  the  right  hand  side  is  the  value  of  the  integral  for  an  ellipsoid.  However,  this 
integral  is  just  equal  to  t*J0.  Thus,  for  types  I,  II,  and  III: 

c  =  ta .  (91) 

Figure  3  shows  the  replacement  ellipsoids  for  two  of  the  appendages  shown  in  Figure  2. 
Figure  1  shows  the  replacement  ellipsoid  for  a  sail. 

Tb  locate  the  appendages  (ie,  obtain  f )  one  requires  the  inveroe  (which  equals  the  trans¬ 
pose)  of  the  B  matrix,  eq.  13: 


Thus: 

£  =  +  fi  cos  O  -  P  sin  O 

y  =  yx -j-psinfl cost  +  Pcosfl cost  (92) 

S  =  zl  +  ji  sin  fl  sin  t  +  P  cos  fi  sin  t . 
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Figure  3  Replacement  ellipsoids  for  two  of  the  appendages  of  Figure  2. 


Unlike  some  methods5,  the  present  theory  does  not  account  for  asymmetric  end  effects 
on  an  appendage.  Theee  effects  are  due  to  the  Bow  around  the  tip  being  different  from  that 
around  the  root.  The  proper  way  of  accounting  for  this  is  by  analysing  the  ffow  around  the 
appendage  plus  its  image  in  the  hull,  and  calculating  the  forces  only  on  that  portion  of  the 
'extended*  appendage  which  protrudes  into  the  Suid.  However,  this  would  have  required  a 
major  alteration  to  the  present  theory,  which  is  as  simple  as  it  is  only  beCaitte  the  added  mass 
coefficients  can  be  related  to  the  total  energy  in  the  Bow.  If  it  were  necessary  to  calculate  the 
force  on  only  a  portion  of  an  ellipsoid,  then  the  complexity  of  the  appendage  analysis  would 
increase  substantially. 

This  criticism  is  moderated  by  the  present  method's  de-emphasis  of  end  effects,  ^he  inte¬ 
grals  of  equations  86  are  only  concerned  with  2-D  aspects  of  the  Sow  norma)  to  the  appendage. 
This  results  in  a  replacement  ellipsoid  with  a  larger  geometric  aspect  ratio  than  the  original 
appendage.  For  example,  the  replacement  ellipsoid  for  an  appendage  with  a  planform 

has  bja  =  1.2. 


9  Optimlxlng  Hull  Ellipsoid  Geometry  for  Interference  Velocity  Calculations 

The  interference  velocity  at  an  appendage  is  determined  by  the  geometry  of  the  hull  replace¬ 
ment  ellipsoid.  In  the  previous  section,  this  geometry  was  chosen  as  the  best  fit  to  the  overall 
hut  geometry.  However,  the  interference  velocity  is  largely  dependent  on  the  local  hull  geom¬ 
etry  opposite  the  appendage.  This  section  describes  the  way  in  which  the  hull  replacement 
ellipsoid  geometry  is  modified  to  beet  represent  the  local  geometry  opposite  the  centroid  of  an 
eppeadage'a  replacement  sllipeoid  prior  to  calculating  the  interference  velocity  at  that  centroid. 

A.  result  of  slender  body  theory  is  that  axial  disturbance  velocities  over  a  slender  hull  are  of 
a  smaller  order  of  magnitude  than  transverse  disturbance  velocities.  The  transverse  velocities 
sr«  primarily  determined  by  the  local  cross-sectional  hull  profile,  while  the  axial  velocities  are 
related  to  the  local  slopes  of  the  longitudinal  profiles.  Ideally,  the  optimised  hull  ellipsoid 
would  model  all  of  these  local  characteristics,  but  compromises  must  be  made.  Since  axial 
interference  velocities  are  email,  no  attempt  is  made  to  match  the  local  longitudinal  profile 
elopes;  such  matching  also  often  results  in  a  dramatic  reduction  in  the  length  of  the  original 
hull  replacement  ellipsoid,  so  that  end  effects  may  then  adversely  affect  the  more  important 
transverse  velocity  predictions.  The  usual  compromise  that  the  croes-eectional  profile  of  the 
hull  ellipsoid  will  be  elliptical,  regardless  of  the  true  hull  profile,  is  also  made. 

The  optimisation  procedure  proceed*  as  follows.  To  improve  the  prediction  of  local  trans¬ 
verse  velocities,  the  transverse  dimensions  of  the  original  hull  replacement  ellipsoid  are  adjusted 
#o  that  the  local  breadth  and  sheer  profile  coordinates  of  the  hull  (obtained  by  linear  interpo¬ 
lation  between  the  specified  profile  coordinates)  are  exactly  matched.  In  general,  this  requires 
adjustment*  to  the  hull  ellipsoid  parameters  6,  c,  and  f .  For  the  longitudinal  direction,  a  is 
only  adjusted  in  those  instances  when  the  appendage  is  too  close  to  the  stern  of  the  hull 
ellipeoid,  as  described  below;  «  is  never  changed. 

In  some  cee*e,  the  aft  end  of  the  original  hull  replacement  ellipeoid  may  be  too  far  forward 
to  allow  the  local  bull  geometry  opposite  the  tail  appendages  to  be  satisfactorily  represented 
with  the  above  modification.  Where  this  occurs,  a  is  usually  increased  in  order  to  make  the 
modified  hull  ellipsoid  terminate  aft  of  the  appendage.  Then  the  transverse  dimensions  of 
ths  sllipeoid  ere  adjusted  to  match  the  hull  profile  coordinates  opposite  the  appendage  (see 
Figure  5  in  the  next  section).  This  usually  provides  a  good  local  fit  at  the  tail  but  results  in  a 
somewhat  narrower  overall  ellipsoid,  which  is  not  a  problem  since  this  modified  hull  ellipsoid 
is  only  used  to  calculate  interference  velocities  at  the  tail. 

If  there  are  cases  where  an  aft  appendage  is  very  close  to  or  even  aft  of  the  end  of  the 
actual  hull,  a  is  only  lengthened  until  the  end  of  the  bull  ellipeoid  is  within  a  distance  d  of  the 
actual  hull  end;  a  is  left  unchanged  if  it  is  already  within  this  distance.  TVial  and  error  has 
shown  that  a  good  value  for  d  is  2  percent  of  the  hull  length.  The  newly  lengthened  (or  not) 
hull  ellipsoid  is  then  matched  to  the  local  hull  profile  coordinates  at  least  a  distance  d  ahead 
of  its  end,  further  forward  if  the  appendage  is  forward  of  this  location. 

If  the  end  of  the  original  hull  replacement  ellipsoid  is  aft  of  the  end  of  the  actual  hull  (a 
very  unusual  situation),  and  if  the  appendage  is  very  close  to  or  aft  of  the  end  of  the  actual 
hull,  no  change  i*  made  to  the  original  hull  ellipsoid. 

At  the  bow,  the  original  hull  replacement  ellipsoid  will  usually  extend  out  beyond  the 
actual  hull,  as  Figure  1  shows.  As  the  presence  of  appendages  very  close  to  the  bow  is  unusual, 
no  attempt  has  been  made  to  optimise  the  hull  ellipsoid  for  such  a  case.  If  an  appendage  is 
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within  a  distance  d  of  or  off  the  forward  end  of  either  the  actual  hull  or  the  original  hull  replace¬ 
ment  ellipsoid,  no  modifications  are  made;  otherwise,  the  transverse  hull  ellipsoid  geometry  is 
optimized  as  for  any  appendage  well  within  the  ends. 

10  An  Example 

A  computer  program  (called  ESAM:  Estimate  Submarine  Added  Meases)  has  been  developed 
which  applies  the  calculations  described  in  this  report  to  submarine  configurations.  Appendix  B 
describes  its  use  and  gives  example  input  and  output  data  files  for  the  generic  submarine 
configuration  shown  in  Figure  4.  Input  data  must  contain  the  hull  sheer  and  breadth  plan 
coordinates,  appendage  geometry,  and  basic  information  about  reference  axes.  The  program 
automatically  calculates  replacement  ellipsoid  geometry  and  outputs  this  information  together 
with  a  complete  set  of  added  mass  predictions.  For  a  complete  submarine  configuration, 
consisting  of  a  hull  and  7  appendages,  the  program  takes  about  2  seconds  of  CPU  time  on  a 
VAX  11/750. 

As  well  as  sketching  the  submarine’s  longitudinal  profiles,  Figure  4  shows  the  size  and 
location  of  the  best  fit  hull  replacement  ellipsoid.  Figure  5  shows  how  this  ellipsoid  is  modified 
by  ESAM  for  the  purpose  of  calculating  interference  velocities  at  the  rudder.  These  replacement 
ellipsoid  geometries,  along  with  those  of  all  the  other  components,  are  listed  in  the  output  file 
of  Appendix  B. 

The  e  values  for  this  hullform  are: 

e,  =  -0.013 ,  es  =  0.013 ,  c4  =  -0.35 ,  ex  =  0.005 
and  for  the  entire  configuration: 

Cy  — -  —0.012  ,  —  0.013 , 

Cjf  —  0.022 ,  €j ^  — -  —0.012  ,  —  0.012  . 

Values  of  e  as  large  as  0.04  were  obtained  when  the  assumption  B{x)/H(x)  =  b/c  was  used 
instead  of  equations  81.  The  asymmetry  in  this  hull  is  more  evenly  distributed  over  its  length, 
and  this  results  in  generally  lower  e  values  than  for  the  previous  hullform;  however,  the  asym¬ 
metry  is  also  quite  large,  so  that  e4  has  a  relatively  large  affect  on  eK . 

Table  1  summarizes  the  added  masses  for  the  complete  submarine.  For  comparison  pur¬ 
poses,  the  table  also  shows  the  predictions  when  hull  interference  velocities  are  not  considered 
and  appendage  thicknesses  are  set  to  zero.  These  conditions  are  all  easily  modelled  by  mak¬ 
ing  simple  parameter  changes  in  the  ESAM  input  file  (interference  effects  are  eliminated  by 
assigning  all  appendages  neutral  orders).  The  prime  over  a  coefficient  indicates  nondimension- 
alization. 

The  following  discussion  deals  with  the  added  mass  terms  on  an  individual  basis.  Note  that 
alternating  coefficients  in  Table  1  must  be  zero  because  of  the  submarine’s  symmetry  about  a 
longitudinal  vertical  plane,  and  since  j/huU  =  0 .  ESAM  does  not  assume  such  symmetry  exists. 

For  this  submarine,  over  20  percent  of  the  magnitude  of  is  due  to  appendage  thickness. 
The  sail  is  most  responsible  for  this;  although  its  section  profile  is  rather  thick  (30  percent), 
its  span  is  not  large  compared  to  some  configurations. 
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a)  Breadth  plan. 


b)  Sheer  plan. 


Figure  4  Generic  submarine  and  hull  replacement  ellipsoid. 


Figure  5  Elevation  view  of  original  and  optimized  hull  replacement  ellipsoids  opposite  the 
rudder;  the  optimized  ellipsoid  is  only  used  for  calculating  interference  velocites 
at  the  centroids  of  the  rudder  replacement  ellipsoids. 


Condition 


k-,m; 


K 


1 

2 

3 

4 

Complete 

Theory 

No 

Interference 

No  Appendage 
Thickness 

Both 

(2)  .  nd  (3) 

-0.000305 

-0.000303 

-0.000236 

-0.000236 

0.0 

0.0 

0.0 

0.0 

-0.00000331 

0.00000278 

-0.00000138 

0.0 

0.0 

0.0 

0.0 

0.0 

0.00000504 

0.00000539 

0.00000087 

0.00000062 

0.0 

0.0 

0.0 

0.0 

-0.0125 

-0.0120 

-0.0126 

-0.0121 

0.0 

0.0 

0.0 

0.0 

-0.0000929 

-0.0000753 

-0.0000989 

-0.0000796 

0.0 

0.0 

0.0 

0.0 

-0.0000805 

-0.0000314 

-0.0000935 

-0.0000395 

-0.00907 

-0.00910 

-0.00904 

-0.00901 

0.0 

0.0 

0.0 

0.0 

-0.0001021 

-0.0000583 

-0.0001130 

-0.0000764 

0.0 

0.0 

0.0 

0.0 

-0.00000463 

-0.00000451 

-0.00000507 

-0.00000493 

0.0 

0.0 

0.0 

0.0 

-0.00000917 

-0.00000680 

-0.00001004 

-0.00000744 

-0.000434 

-0.000423 

-0.000435 

-0.000422 

0.0 

0.0 

0.0 

0.0 

-0.000541 

-0.000525 

-0.000543 

-0.000526 

Table  1  Added  mass  predictions  for  the  submarine  of  Figure  4.  Column  1  presents  the 
predictions  of  the  complete  theory  (component  contributions  to  these  values  are 
shown  in  the  output  Sle  of  Appendix  B);  column  2  neglects  interference  effects; 
column  3  sets  all  appendage  thicknesses  to  zero;  column  4  neglects  both  interference 
effects  and  appendage  thicknesses.  Coefficients  are  with  respect  to  an  origin  located 
on  the  bull  centerline  opposite  the  centroid  of  the  hull  replacement  ellipsoid. 
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Table  1  shows  how  X ^  is  entirely  dependent  on  the  modelling  of  both  interference  effects 
and  appendage  thickness.  The  present  theory  models  the  hull  with  an  ellipsoid,  which  has 
horizontal  and  vertical-transverse  planes  of  symmetry  not  matched  by  the  actual  hull.  It 
predicts  that  the  hull  makes  no  contribution  to  X ^  when,  in  fact,  one  would  expect  a  small 
contribution  to  be  present.  The  size  of  such  a  contribution  is  unknown.  However,  the  sensitivity 
of  X^  to  small  asymmetries,  such  as  the  negative  sweepback  of  the  sail  (Qifcil  =  -3.3°),  can 
be  determined:  when  ng,sl  is  set  to  zero,  the  magnitude  of  changes  by  57  percent.  For  this 
submarine  configuration,  then,  the  current  X^  prediction  may  only  be  an  order  of  magnitude 
estimate.  For  configurations  with  hullforms  more  closely  exhibiting  the  symmetries  of  an 
ellipsoid,  the  prediction  of  X^  might  be  better. 

Xj  is  essentially  determined  by  sail  thickness,  with  interference  effects  being  of  secondary 
importance.  The  contribution  to  X^  from  the  hull  (which  is  nonzero  because  the  hull  centroid 
is  displaced  upw  Is  from  the  hull  centerline)  is  an  order  of  magnitude  less  than  that  from  the 
sail. 

Note  that  4  percent  of  the  value  of  Y i  is  due  to  interference  effects.  These  effects  increase 
the  sail’s  contribution  to  Yi  by  85  percent.  Table  2  compares  the  present  theory’s  predictions 
of  the  sail’s  contribution  to  Y ^  to  the  theories  mentioned  in  the  Introduction  (which  account 
for  neither  hull  interference  velocities  nor  appendage  thickness;  Aucher  models  the  sail’s  image 
in  the  hull  while  Humphreys  and  Watkinson  do  not). 


Present  Theory 

Humphreys 

complete 

without 

neither  thickness 

and 

Aucher5 

interference 

nor  interference 

Watkinson4 

9.7 

5.2 

5.8 

5.2 

6.0 

Table  2  The  contribution  to  104Y0'  of  the  B&il  of  Figure  4. 


YffNfi  and  Z^, are  difficult  to  predict  accurately  if  they  are  referenced  to  a  point 
close  to  the  center  of  buoyancy  (CB).  This  is  because  they  are  being  estimated  iu  the  region  of 
their  zeros  (which,  incidently,  are  only  approximately  at  the  CB).  Of  more  importance  is  the 
prediction  of  Y^/Y^ ,  Ni)/Yi ,  -Z^Z^ ,  and  -M^/Z^ ,  the  longitudinal  locations  of  the  zeros 
of  the  coefficients. 

Although  Table  1  shows  very  little  change  in  Z ^  for  the  various  conditions,  somewhat 
larger,  partially  compensating  changes  are  occuring  in  the  component  contributions  to  it.  Fo- 
example,  properly  accounting  for  interference  effects  (ie,  going  from  column  2  to  column  1) 
roughly  cuts  in  half  the  contributions  from  the  sail  (due  to  its  thickness)  and  the  sailplanes, 
while  at  the  same  time  increasing  by  a  third  the  contributions  from  the  sternplanes.  The  effects 
of  accounting  for  thickness  and  interference  are  more  easily  appreciated  in  the  changes  that 
result  in  M 4  ,  since  is  sensitive  to  the  longitudinal  location  at  which  these  changes  take 
place. 

The  K  derivatives  are  dominated  by  the  presence  of  the  sail.  Although  the  hull  accounts 
for  a  third  of  the  value  of  K 0  ,  it  accounts  for  less  than  10  percent  of  Kp  and  makes  no  con¬ 
tribution  to  K+.  The  cross-flow  velocities  associated  with  v  and  f  result  in  large  interference 
velocity  effects  in  the  seal’s  contribution  to  and  K+ ,  but  not  as  large  as  one  might  expect 
(cf,  equations  29). 


11  Conclusions 


A  theory  for  predicting  all  the  added  mass  coefficients  for  a  multi-component,  deeply  submerged 
vehicle  has  been  presented.  A  program  applying  this  theory  to  conventional  submarine  shapes 
gives  very  rapid  numerical  predictions  of  the  coefficients.  In  keeping  with  the  way  most  current 
experimental  investigations  of  these  characteristics  are  carried  out,  the  theory  ignores  both  the 
incidence  of  the  vehicle  to  the  oncoming  flow  and  time  in  an  unsteady  flow. 

Added  masses  we  calculated  by  summing  up  the  contributions  from  replacement  ellipsoids 
optimally  sized  and  positioned  to  best  represent  each  body  component.  The  hull  replacement 
ellipsoid  is  given  special  attention;  it  is  sized  in  such  a  way  that  errors  associated  with  the 
replacement  we  minimized. 

The  theory  accounts  for  both  appendage  thickness  and  hull  interference  velocities  at  ap¬ 
pendages.  Modelling  these  chwacteristics  is  shown  to  be  necesswy  if  reasonable  predictions 
of  the  ‘off-diagonal’  added  mass  coefficients,  responsible  for  at  least  some  degree  of  coupling 
between  motions  in  different  degrees  of  freedom,  are  to  be  made.  In  addition,  these  chwacter¬ 
istics  we  responsible  for  changes  to  five  of  the  principal  (ie,  ‘diagonal’)  added  mass  coefficients 
of  from  1  to  6  percent;  modelling  appendage  thickness  may  increase  the  prediction  of  X u  by 
as  much  as  30  percent. 

If  more  accuracy  than  the  present  theory  provides  is  required,  the  theory  should  progress 
on  at  least  two  fronts  simultaneously:  circulation  should  be  properly  accounted  for  in  the  flow, 
so  that  the  added  masses  become  a  function  of  the  state  vwiables;  and,  a  better  geometri¬ 
cal  representation  of  the  multi-component  vehicle  should  be  developed.  Since  added  mass  is 
merely  a  measure  of  the  kinetic  energy  in  the  flow,  properly  accounting  for  boundary  layer 
growth,  sepwation,  and  the  flow  structures  they  generate,  all  of  which  effectively  use  up  or 
re-distribute  kinetic  energy,  may  also  be  necessary.  The  effect  of  the  propeller  may  also  need 
to  be  considered.  Accounting  for  any  of  these  features  would  require  a  major  extension  of  the 
present  theory. 

Future  work  should  consider  the  effect  of  a  free  surface  on  the  added  masses.  As  well, 
sensitivity  studies  should  be  cwried  out  in  order  to  establish  the  accuracy  with  which  the 
added  mass  coefficients  should  be  known. 
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APPENDIX  A 


Elliptic  Functions 


This  appendix  describes  methods  for  evaluating  elliptic  integrals  of  the  first  and  second  kinds. 
Extensive  use  is  made  of  the  theory  of  theta  functions  and  Jacobian  elliptic  functions,  as 
presented  by  Whittaker  and  Watson11,  Chapters  21  and  22,  and  Watt18,  Appendix  B.  The 
serious  reader  is  advised  to  at  least  look  over  Reference  11. 

The  fundamental  elliptic  integral  of  the  first  kind  is: 


[* _ dt _ 

U~Jo  y/{l  -  t*)(l  -  kH*) 


=/, 


de 


o  y/l  -  k*  sin2  0 


(Al) 


where  k  is  the  modulus  of  the  integral:  0  <  k  <  1.  These  integrals  are  related  through  the 
variable  transformation  t  =  sin  0 .  If  z  =  1 ,  so  that  <j>=  jt/2,  the  integral  is  called  a  complete 
elliptic  integral  of  the  first  kind  and  is  represented  by  the  symbol  K .  The  complementary  mod* 
ulus,  k'  =  \/l  -  fc3,  determines  the  only  other  complete  elliptic  integral  of  the  first  kind,  K' , 
which  is  the  same  function  of  k'  as  K  is  of  k . 


The  Jacobian  elliptic  functions  sn(u,  k) ,  cn(u,  k) ,  and  dn(u,  k)  (usually  written  simply  as 
snu,  cnu,  and  dn u  if  the  value  of  the  modulus  to  be  used  is  k)  are  related  to  each  other  via: 


sn3  u  +  cn3  u  =  1 ,  dn3  u  +  fc3  sn2  u  =  1 . 


(A2) 


Also: 


—  sn  u  =  cn  u  dn  u  , 
au 


d  d  , 

—  cnu  =  -snudnu,  —  dnu  =  -Arsnucnu. 
au  au 


Substituting  t  =  sn  u  in  the  integrand  of  equation  Al  gives: 

u  =  sn-1  x  =>  x  —  snu 


(A3) 


(A4) 


so  that  the  inverse  sn  function  can  be  defined  as  an  incomplete  elliptic  integral  of  the  first  kind. 
Whereas  the  trigonometric  functions  sin  u  and  cos  u  are  singly  periodic  in  the  complex  u  plane 
(they  are  periodic  along  the  real  axis  but  not  along  the  imaginary  axis),  the  Jacobian  elliptic 
functions  are  doubly  periodic.  For  example,  snu  =  sn(u ± 4mA ± 2niK')  for  m,n  =  0, 1,2, .. ., 
leading  to  cells  of  periodicity  in  the  complex  u  plane. 

The  fundamental  elliptic  integral  of  the  second  kind  is: 


dn3(u,fc)  du 


(A5) 


or,  making  the  substitution  t  =  sn  u ,  followed  by  t  =  sin  9 : 


E(u)=  r^=  J’j'dt  =  Vl  —  k*  sin3  6  dO . 


(A6) 


When  *  =  1,  the  integral  is  designated  E  (ie,  E  =  E(K))  and  is  called  a  complete  elliptic 
integral  of  the  second  kind.  E '  is  the  only  other  complete  elliptic  integral  of  the  second  kind 
and  it  is  the  same  function  of  k'  as  E  is  of  k . 

When  numerical  work  must  be  done,  the  use  of  theta  functions  is  advised.  The  Jacobian 
elliptic  functions  can  be  constructed  from  the  four  quaai-doubly  periodic  theta  functions: 


sn(u,  k) 


1  *i(*>g) 

y/k#A(*,q)  ’ 


cn(u,  k) 


v  **4(*,?r 


dn(u,  Jfc)  =  v/*7 


where  the  theta  functions  are  more  naturally  defined  in  terms  of: 


T 


(A7) 


(A8) 


The  parameter  q  is  purely  a  function  of  k ,  and  is  to  the  theta  functions  what  k  is  to  the 
Jacobian  elliptic  functions: 

q  =  e-*K'/K.  (A9) 

Whittaker  and  Watson  give: 


t?l(*)  =  2 q*  sin  *  -  2 q*  sin  3*  +  2 q3^  sin5z  -  •  •  • 
i?a(z)  =  2 q*  cos*  +  2g«  cos 3*  +  2 q3*  cos 5*  -+-••• 
t?3(*)  =  1  -r  2 q cos  2*  +  2 q4  cos 4*  +  2 q°  cos 6*  +  •  •  • 
t?4(z)  =  1  —  2qcos2z  +  2q4  cos  4*  -  2g*cos6z  +  •  •  •  . 


(A10) 


Now,  if: 


1  -  \[k< 

i  +  v/ifc7  (i  +  kr)  (i  •+•  v*T 


(All) 


(the  latter  form  of  this  equation  being  the  most  accurate  for  numerical  calculations  when  k  is 
very  small)  then  Whittaker  and  Watson  show  that: 


q  =  e0  +  2e£  +  15e»  +  150cJ9  +  0(l000ej7). 


(A12) 


</,c o  0  as  A:  -♦  0  and  q,2e0  — ►  1  as  k  — »  1.  Nevertheless,  this  equation  is  accurate  for 
surprisingly  large  values  of  k .  For  complete  coverage  of  the  range  of  k  values  right  through  to  1, 
however,  one  replaces  k1  with  k  in  equation  All  and  then  calculates  q1  from  equation  A12. 
q  and  q'  are  related  via: 

?/  =  e-^//C'  =eKa/ln,  (A13) 

With  q  determined,  equations  A14  and  A15  give  the  complete  elliptic  integrals  of  the  first  and 
second  kinds: 

—  (l  4-  2q  +  2  q*  +  2  q9  +  2q 16  +  •  •  •) 

?-?(?)-  m 

Legendre’s  relation  gives  E': 


EK'+E'K  -  KK' 


(A15) 


I 


where  use  has  been  made  of  a  manifestation  of  Jacobi's  imaginary  transformation: 


'  cn(u ,k) 


>(*«,  *')  =  — p-rr ,  dn(tu,  k')  =  dP|U>  Q 
cn («,«)  cn(u,k) 


Then,  with: 


it  follows  that,  as  q'  —>  0  : 


2F,U’ 


(A21) 


(A22) 


cosh  (dF2u)  ~ f/ 11  “  9 ~  4f,S)  +  ?,*(3  “  20f,J  +  32f,<) 

-  g'lf(8  -  76f'*  +  272f'4  -  320f'®)  +  0(l000f,8g'18)] . 


(A23) 


In  the  range  0  <  u  <  K,  f'  varies  from  a  minimum  of  approximately  1  at  u  =  0  to  a 
maximum  of  1/2 q'  at  u  =  K .  Thus,  equation  A23  will  not  be  as  accurate  as  its  counterpart, 
equation  A19;  however,  between  the  two  of  them,  the  complete  range  of  k  values  can  be 
accurately  covered. 

The  Incomplete  Elliptic  Integral  of  the  Second  Kind 

Whittaker  and  Watson  show  that  the  incomplete  elliptic  integral  of  the  second  kind  can  be 
written: 

E{u)  ~  u  ^  j K  (A24^ 

where  t?4(z)  is  the  derivative  of  i?4(z)  with  respect  to  z .  This  expression  can  now  be  evaluated 
using  equations  A10.  However,  there  is  a  more  efficient  and  accurate  method  available. 

By  taking  the  derivative  of  the  logarithm  of  the  first  of  equations  A16  and  the  last  of 
equations  A7,  Watt  shows  that: 


d4(z)  _  2 K_  fc*snucnu  t?3(2z,g4) 

04(*)  *  dnu  +  >A7+  »M2*>?4) 


(A25) 


Thus,  if  the  Jacobian  elliptic  functions  of  u  are  known,  the  number  of  terms  that  need  to  be 
evaluated  in  a  calculation  of  E( u)  can  be  substantially  reduced.  This  is  particularly  true  if  u 
is  imaginary,  as  in  the  next  paragraph. 

Equation  A25  works  well  if  q  <£  1 .  If  it  iB  not,  another  form  of  Jacobi’s  imaginary 
transformation  can  be  exploited: 


.  snudnu  ... 

Mu)  =  u  4- - -I-  tE( »u, k1) . 

'  cnu 

E(iu,k')  is  then  expanded  in  a  manner  analogous  to  equation  A24: 


E(iu  Jfc'l  -  iu—  +  * 

E{xu,k)~  tuK,  +  2K'd<(z',q')' 


(A26) 


(A27) 


Hence,  for  q1  1 : 


PM  ufl  E>\  l  TO“(*,cn«  +  V'*dnu)  .  ..  «■  «?s(2*',g'4)  f.„av 

z;(«)  _  „  +  _____  +  .  (AJ8) 
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Program  ESAM  User  Instructions 


This  »ppendix  gives  the  user  instructions  and  describes  the  data  input  and  output  formats  for 
ESAM  (Estimate  Submarine  Added  Masses),  a  FORTRAN  77  program  which  performs  the 
calculations  required  by  the  theory  developed  in  the  main  body  of  this  report.  After  compiling, 
linking,  and  run  ing  the  program,  the  user  is  asked  to  give  the  name  of  the  input  file  which 
must  contain  all  the  data  for  a  complete  run.  The  user  simply  enters  this  file  name;  except  for 
possible  error  messages  being  written  to  the  designated  output  unit  (the  terminal  cn  a  VAX), 
there  is  no  other  interaction  with  the  user.  Output  is  sent  to  a  newly  created  file,  ESAM.OUT. 
If  for  some  reason  the  program  terminates,  ESAM.OUT  vill  show  how  many  components  were 
successfully  processed  before  termination  occurred.  If  an  error  message  is  written  to  the  screen 
by  ESAM,  it  will  also  be  written  to  ESAM.OUT. 


B.l  Input  Format 

This  section  describes  the  data  (and  the  order  in  which  they  must  occur)  that  the  input  file 
must  contain.  In  general,  data  should  appear  exactly  as  specified  in  the  input  format  summary; 
spaces  should  not  be  left  between  lines.  However,  since  READ  statements  in  ESAM  are  free 
format  (except  for  those  reading  character  constants),  numerical  data  can  be  split  between 
different  lines  if  desired.  ESAM  does  its  best  to  test  the  input  data  to  make  rvre  the  rules 
described  below  are  being  followed;  if  they  are  not,  an  appropriate  error  message  is  given. 

Input  format  summary 


NC 

XDIR 

XOFF 

ZOFF 

ELL 

TYPE  COMMENT 

XSHIFT ,  YSHIFT ,  ZSHIFT  } 


XO.B(XO) ,H(XO) .HBAR(XO)  ' 
XI ,B(X1) ,H(X1) ,HBAR(X1) 

XN.B(XN) ,H(XN) .HBAR(XN) 


Only  need  to 
enter  XI  .B(Xi) 
if  profile  is 
designated 
axiSymmetric. 


> 


TYPE  COMMENT  1 
XI ,Y1 ,Z1 ,PHI .XH2.YH2 .XH3.YH3 ,XH4 ,YH4 ,T  J 
TYPE  COMMENT  1 
XI .Yl.Zl pPHI,Cl.XH4,YH4,C4,XH6.YH6,C5.T J 
TYPE  COMMENT  ) 
A.B.C.XBAR.YBAR.ZBAR.PHI.OMEGA  j 


assuming  type  is  H. 


assuming  type  is  1  or  2. 
assuming  type  is  3. 
assuming  type  is  E. 


Various  TYPEs  of  components  are  added  to  this  list,  as  required. 


Definition  of  variable* 

NC  An  integer  specifying  the  number  of  components  which  make  up  the  submarine  and 

for  which  data  is  about  to  be  specified.  Array  dimensioning  limits  NC  to,  at  most,  20. 

XDIR  The  added  mass  coefficients  output  from  this  program  are  a /ways  calculated  rela¬ 
tive  to  the  conventional  body  fixed  system  of  axee  used  by  the  DTNSRDC  standard 
submarine  equations  of  motion*.  Although,  in  general,  this  will  also  be  the  coor¬ 
dinate  system  used  to  input  submarine  geometry,  provision  is  made  for  input  data 
x  coordinates  that  increase  opposite  to  the  standard  sense.  Note  that  all  data  input 
during  a  run  of  the  program  must  be  referenced  to  an  identical  set  of  input  axes, 
whatever  they  may  be. 

XDIR  is  a  character  specifying  the  direction  of  increasing  x  for  all  input  data:  if 
XDIR  is  P  (for  positive),  the  coordinate  system  for  the  input  data  will  be  the  standard 
body  fixed  system  of  axes,  with  x  coordinates  increasing  from  stern  to  bow;  if  XDIR 
is  N  (for  negative),  the  input  coordinate  system  will  still  be  the  standard  body  axes 
except  that  the  x  direction  will  be  reversed  so  that  x  coordinates  increase  from 
bow  to  stern  (y  and  n  coordinates  always  increase  going  to  starboard  and  going 
downwards,  respectively,  regardless  of  XDIR  being  P  or  N). 

The  summary  of  geometrical  information  written  to  ESAM.OUT  is  also  con¬ 
verted  to  the  standard  axes  used  by  the  added  mass  coefficients.  The  next  paragraph 
specifies  how  the  origin  of  this  system  is  determined. 

XOFF  All  data  must  be  input  relative  to  the  same  origin,  which  can  be  located  anywhere 

ZOFF  that  is  convenient.  XOFF  and  ZOFF  (x-offset  and  z-offset)  are  real  constants  that  are 

subtracted  from  the  input  x  and  z  coordinates  to  give  a  new  origin  for  the  body 
fixed  axes  to  which  all  output  is  referenced.  Note  that  the  input  and  output  data 
have  identical  y  axis  origins. 

Sometimes  it  is  desirable  to  locate  the  output  axes  origin  at  the  center  of 
bouyancy,  even  though  knowledge  of  its  location  is  unknown  beforehand.  Although 
ESAM  does  not  accurately  calculate  this  point,  the  location  of  the  centroid  of  the 
replacement  hull  ellipsoid  closely  approximates  it.  If  the  user  wishes  to  use  either 
or  both  of  the  x  and  z  coordinates  of  this  location  for  specifying  his  output  origin, 
he  need  only  input  C  (for  centroid)  for  either  or  both  of  XOFF  and  ZOFF  If  this  is 
done,  the  centroid  of  the  first  component  entered  locates  the  new  or1' gin,  so  this  first 
component  must  represent  the  hull. 

ELL  A  real  constant  equal  to  the  hull  length;  its  magnitude  is  used  to  nondimensionalize 
the  added  mass  coefficients.  If  ELL  =  0,  nondimensionalisation  does  not  take  place; 
in  this  case,  the  added  masses  need  only  be  multiplied  by  the  fluid  density  (in  units 
consistent  with  those  of  the  input  data)  to  obtain  the  actual  added  masses. 

ELL  is  also  used  to  determine  the  format  in  which  the  added  mass  coefficients  are 
written  to  ESAM.OUT:  if  ELL  is  greater  than  zero,  the  coefficients  are  written  via 
the  edit  descriptor  F15 . 10;  if  ELL  is  less  than  or  equal  to  zero,  the  descriptor  1PE15 . 8 
is  used.  Thus,  the  user  chooses  either  F  or  E  formatted  output  for  nondimenaionalized 
coefficients  by  making  ELL  positive  or  negative.  Dimensional  coefficients  are  always 
output  in  E  format. 
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Subsequent  records  contain  component  information.  The  rules  given  below  will  determine  the 
order  in  which  some  components  should  be  entered. 

TYPE  This  is  a  character  constant  of  length  2  and  it  precedes  the  input  format  for  all 
TTPEs  of  components.  The  first  character  specifies  the  component  type.  Component 

types  are: 

H  (for  Hull), 

1  (for  type  I  appendage), 

2  (for  type  II  appendage), 

S  (for  type  m  appendage),  and 

E  (for  Ellipsoid). 

The  second  character  specifying  the  TYPE  is  the  order.  It  is  used  to  determine  if 
interference  effects  are  to  be  accounted  for  on  this  component.  There  are  3  orders: 

P  (for  Primary), 

S  (for  Secondary),  and 
N  (for  Neutral). 

A  primary  component  can  only  be  of  type  H  or  E;  it  determines  the  interference 
velocity  field  at  all  secondary  components  entered  between  it  and  the  next  primary 
component;  its  added  masses  are  not  corrected  for  interference  effects  from  any 
other  component.  Although  ESAM  does  not  prevent  the  user  from  entering  more 
than  one  primary  component,  to  do  so  makes  little  sense  since  it  implies  more  than 
one  independent  configuration  is  being  processed  at  once,  and  ESAM  sums  all  the 
component  added  masses  into  one  total  added  mass. 

A  secondary  component  can  be  associated  with  any  type;  it  must  be  preceded 
by  a  primary  component. 

A  neutral  component  can  also  bo  associated  with  any  type;  it  has  no  effect  on 
any  other  component;  it  is  not  corrected  for  interference  effects;  it  can  be  entered 
anywhere  amongst  the  other  components. 

There  can  be  any  number  of  primary,  secondary,  and  neutral  order  components 
during  a  run,  up  to  a  total  of  20  (limited  only  by  array  dimensioning).  Of  course, 
primary  and  secondary  components  must  be  ordered  in  a  meaningful  way. 

Examples  of  various  TYPEs  are:  HP,  HS,  HN,  IS,  SN,  and  EP.  The  first  8  compo¬ 
nent  designations,  here,  indicate  hull  components  which  are  of  primary,  secondary, 
and  neutral  orders,  respectively.  The  fourth  component  is  a  type  1  appendage  of 
secondary  order.  The  fifth  is  a  type  3  appendage  of  neutral  order.  The  last  one 
specifies  an  ellipsoid  of  primary  order. 

COMMENT  A  character  constant  of  length  12;  it  must  be  on  the  same  line  as  TYPE  and 
separated  from  it  by  2  characters.  COMMENT  is  any  phrase  the  user  wants  to  use 
to  describe  the  associated  component;  it  is  only  used  for  descriptive  purposes  in 
ESAM.OUT,  and  may  be  left  blank.  Actually,  ESAM  reads  in  ‘TYPE  COMMENT’  as  a 
single  constant  of  length  16. 
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If  the  TYPE  designation  apecifiee  type  H,  then  ♦  end  ft  are  set  to  aero  and  the  following  input 

format  mutt  follow  TYPE,  regardleee  of  the  order  : 

XSHIFT.YSHIFT.ZSHIFT  Real  conitanta  specifying,  in  the  usual  input  coordinate  system, 

a  shift  to  be  applied  to  the  hull  component  coordinates.  XSHIFT 
and  Z8HIFT  are  effectively  added  to  the  values  of  Xi  and  HBAR(Xi)  defined  below. 
9  is  equated  to  Y8HXFT.  Although  theee  data  are  unnecessary  for  a  primary  hull 
component,  which  can  make  use  of  the  XOFF  and  ZQFF  inputs  and  for  which  f)  is 
invariably  sero,  they  are  convenient  for  secondary  components  which  are  not  hulls 
but  instead  are  pods  or  some  other  cylindrical  type  of  body  mounted  out  from 
the  hull.  If  the  shift  coordinates  are  used  with  a  primary  component,  one  should 
ensure  that  subsequent  secondary  components  are  properly  positioned  relative  to 
the  shifted  primary  component;  unlike  XOFF  and  ZOFF,  the  shift  coordinates  are 
only  applied  to  the  hull  component  they  are  associated  with. 

Since  use  of  the  shift  coordinates  will  usually  be  unnecessary,  this  record  may  be 
bypassed,  in  the  usual  way  for  list  directed  input,  by  replacing  it  with  a  /.  Bypassing 
the  record  is  equivalent  to  setting  the  shift  coordinates  to  sero. 

N  An  integer  specifying  the  number  of  sections  the  component  is  divided  into  for  the 

purposes  of  numerically  specifying  its  profile.  Note  that  N  sections  are  delineated 
by  N  +  1  stations  and  that  profile  coordinates  must  be  given  at  each  station.  Array 
dimensioning  restricts  N  to  being  less  than  or  equal  to  100. 

XO.B(XO)  ,H(XO)  .HBAR(XO)  Here,  Xi  is  the  x  coordinate  of  the  1th  axial  station;  it  must 

XI  .8(11)  ,H(X1)  .HBAR(Xl)  not  decrease  as  its  index  increases.  Xi  will  increase  from  stern 

to  bow  if  XDIR  is  positive  and  from  bow  to  stern  if  XDIR  is 

XN.B(XN)  ,H(XN)  .HBAR(XN) 

B(Xi)  is  the  hull  component’s  local  breadth;  it  is  assumed 
to  be  symmetrical  about  y  =  g . 

H(Xi)  is  the  local  height  as  seen  in  a  sheer  plan;  it  is  not  assumed  to  be 
symmetrical  in  any  sense. 

HBAR(Xi)  is  the  x  coordinate  of  the  mid-point  of  the  local  height,  H(Xi). 
HBAR(Xi)  can  be  used  to  describe  two  things:  1)  the  “camber*  of  the  elevation 
profile  and  2)  the  height  of  the  component  above  or  below  the  input  x  axis.  Al¬ 
though  this  latter  capability  allows  the  component  to  be  raised  or  lowered  relative 
to  the  input  x  axis,  the  simplest  way  to  do  this  is  to  use  ZSHIFT. 

Note  that  the  N  +  1  stations  need  not  be  evenly  distributed  over  the  hull  com¬ 
ponent.  Profile  coordinates  can  be  concentrated  in  areas  of  high  curvature,  such  as 
the  nose,  and  eliminated  from  regions  of  sero  curvature  (ie,  where  the  breadth  and 
sheer  profiles  are  straight  lines).  This  is  because  the  only  approximation  made  in  the 
numerical  integration  of  the  integrals  of  equations  74  and  83  is  that  the  input  profile 
coordinates  are  all  joined  by  straight  lines;  if  in  fact  they  are,  then  the  integrals  are 
evaluated  exactly.  Thus,  a  hull  component  with  a  long  cylindrical  mid-body  section, 
for  example,  only  needs  to  have  the  mid-section  profile  specified  at  its  endpoints. 

For  axisymmetric  profiles,  only  Xl,B(Xi)  need  be  entered.  However,  ESAM 
must  first  be  told  that  this  will  be  the  case;  this  is  done  by  appending  an  S  (for 
axiSymmetric)  to  the  TYPE.  For  an  axisymmetric  hull  component  of  primary  order, 
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‘TYPE  COMMENT’  might  look  like  'HP8  HULL  with  only  1  character  sepa¬ 

rating  HPS  and  HULL.  For  the  axieymmetric  case,  ESAM  always  looks  for  XI  ,B(Xi) 
at  the  beginning  of  a  new  record.  Thus,  if  a  user  wants  to  specify  an  axisymmet- 
ric  hullform  using  an  asymmetric  data  file,  deletion  of  the  H(Xi)  ,HBAR(Xi)  data  is 
unnecessary  if  it  is  on  the  same  lines  as  the  Xi,B(Xi)  data.  Designating  the  hull 
component  as  axisymmetric  has  two  other  minor  consequences:  1)  ESAM’s  calcula¬ 
tions  will  be  a  little  quicker,  since  certain  simplifications  can  be  made;  2)  the  errors 
specifying  the  accuracy  with  which  the  replacement  ellipsoid  is  able  to  model  the 
true  hull  component  are  not  written  to  ESAM. OUT,  since  they  are  always  sero  for 
an  axisymmetric  profile. 

It  is  permissible  to  input  axisymmetric  hull  components  using  the  asymmetric 
format. 

This  completes  the  hull  component  input  data.  Now,  the  input  file  will  either  terminate  or  the 
next  entry  will  be  another  TYPE  variable,  followed  by  its  associated  component  information. 

If  the  TYPE  designation  specifies  a  type  1  or  2  appendage,  then  the  following  input  format 
must  follow  TYPE.  The  user  is  reminded  that  the  only  difference  between  a  type  1  and  2 
appendage  is  in  how  the  sweepback  angle,  fl,  is  calculated:  for  type  1,  fi  equals  the  sweepback 
angle  of  the  bisector  of  the  angle  formed  by  the  leading  and  trailing  edges;  for  type  2,  fl  is 
the  angle  between  the  line  joining  the  leading  and  trailing  edge  root  coordinates  and  the  body 
x  axis. 

X1,Y1,Z1.PHIIXH2,YH2,XHS,YHS.XH4IYH4>T  These  are  all  real  constants.  Because  of  the 

free  format  HEAD  statement,  they  can  be 
entered  on  different  lines  if  the  user  wishes. 

XI  ,Y1 .21  are  the  coordinates  of  the  appendage  trailing  edge  at  the  root,  using 
the  usual  input  coordinate  system. 

PHI  is  the  angle  4  (in  degrees)  the  y  axis  must  rotate  about  the  x  axis  to 
bring  it  into  the  plane  of  the  control  surface;  for  a  sail,  for  example,  PHI  could  be 
input  as  -00  or  270. 

XH2.YH2  are  the  coordinates  (xs,yj)  of  the  leading  edge  nt  the  root,  using  an 
unswept  local  coordinate  system  with  its  origin  at  (Xl.Yl.Zl),  the  XH  axis  point¬ 
ing  towards  the  bow  (parallel  to  the  body  x  axis),  the  YH  axis  pointing  outwards 
towards  the  appendage  tip,  and  with  the  ZH  axis  perpendicular  to  the  planform  of 
the  appendage  (see  the  text,  Figure  2). 

XH3.YH3  are  the  local  coordinates  (fs, ys)  of  the  leading  edge  at  the  tip. 

XH4.YH4  are  the  local  coordinates  (£4,^4)  of  the  trailing  edge  at  the  tip. 

T  is  the  thickness  to  chord  ratio,  t ,  of  the  section  profile. 

Note  that  triangle*  can  be  input  in  this  format.  If  (XH3.YH3)  is  coincident 
with  (XH4.YH4),  processing  continues  normally.  If  (XH2.YH2)  is  coincident  with 
(XH3.YH3)  or  if  (XH4.YH4)  is  coincident  with  (0,0),  then  the  appendage  is  automat¬ 
ically  converted  to  a  type  2  appendage,  if  it  isn’t  already,  so  that  fl  is  taken  to  be 
the  argument  of  the  (XH2.YH2)  vector. 

This  completes  the  type  1  and  2  appendage  component  input  data.  Now,  the  input  file  will 
either  terminate  or  the  next  entry  will  be  another  TYPE  variable,  followed  by  its  associated 
component  information. 


If  the  TYPE  designation  specifies  a  type  3  appendage,  then  the  following  input  format  must 
follow  TYPE.  The  user  is  reminded  that  the  sweepback  angle,  Cl ,  is  taken  to  be  zero  for  this 
appendage.  As  noted  in  the  text,  a  pair  of  appendages  separated  by  the  hull  should  not  be 
represented  with  a  type  3  appendage  if  one  wishes  to  account  for  interference  effects. 

XI ,  Y1 ,  Z1 , PHI .  Cl ,  XH4 .  YH4 ,  C4 ,  XH5 ,  YH5 .  C6 .  T  These  are  all  real  constants;  they  can  be 

entered  on  different  lines  if  the  user  wishes. 

XI ,  Y1  ,Z1  are  the  coordinates  of  the  appendage  trailing  edge  at  the  root,  using 
the  usual  input  coordinate  system. 

PHI  is  the  angle  #  (in  degrees)  the  y  axis  must  rotate  about  the  x  axis  to 
bring  it  into  the  plane  of  the  control  surface;  for  a  sail,  for  example,  PHI  could  be 
input  as  -00  or  270. 

Cl  is  the  root  chord  length  c1 ,  assumed  parallel  to  the  body  x  axis. 

XH4.YH4  are  the  coordinates  (x4,y4)  of  the  trailing  edge  at  the  upper  tip, 
using  an  inswept  local  coordinate  system  with  its  origin  at  (XI ,  Y1 ,  Zl),  the  XH  axis 
pointing  towards  the  bow  (parallel  to  the  body  x  axis),  the  YH  axis  pointing  towards 
the  upper  tip,  and  with  the  ZH  axis  perpendicular  to  the  planform  of  the  appendage. 

C4  is  the  upper  tip  chord  length  c4 ,  parallel  to  the  body  x  axis. 

XH5.YH5  are  the  local  coordinates  (xs,y5)  of  the  trailing  edge  at  the  lower  tip 
(so  YH5  will  be  negative). 

C5  is  the  lower  tip  chord  length  c6 ,  parallel  to  the  body  x  axis. 

T  is  the  thickness  to  chord  ratio,  t ,  of  the  section  profile. 

In  this  discussion,  the  terms  ‘upper’  and  ‘lower’  tip  were  used  to  orient  the 
reader.  These  are  really  misnomers,  however,  since  PHI  can  be  used  to  rotate  the 
planform  to  any  desired  angle. 

This  completes  the  type  3  appendage  component  input  data.  Now,  the  input  file  will  either 
terminate  or  the  next  entry  will  be  another  TYPE  variable,  followed  by  its  associated  component 
information. 

If  the  TYPE  designation  specifies  a  type  E  component,  then  the  following  input  format  must 
follow  TYPE.  This  component  is  an  ellipsoid  which  the  user  can  place  anywhere  he  likes.  Since 
it  is  already  an  exact  ellipsoidal  representation  of  itself,  ESAM  does  not  replace  it  with  any 
other  geometry.  If  this  component  is  of  primary  order,  then  PHI  and  OMEGA  must  be  zero. 

A.B.C.XBAR.YBAR.ZBAR.PHI.OMEGA  These  are  all  real  constants;  they  can  be  entered  on 

different  lines  if  the  user  wishes. 

A,B,C  are  the  ellipsoid  principal  semi-axis  lengths  a,b,c:  A  is  associated  with 
the  x'  axis  (the  x  axis  if  OMEGA  =  0),  B  with  the  y'  axis  (the  y  axis  if  PHI  = 
OMEGA  =  0),  and  C  with  the  e'  axis  (the  z  axis  if  PHI  =  OMEGA  =  0). 

XBAR ,  YBAR ,  ZBAR  are  the  coordinates,  using  the  usual  input  data  coordinate 
system,  of  the  ellipsoid’s  centroid. 

PHI ,  OMEGA  are  the  angles  $  and  Cl  (in  degrees)  defining  the  orientation  of  the 
ellipsoid  principal  axes  relative  to  the  body  fixed  axes:  PHI  is  the  angle  the  y;  axis 
has  rolled  away  from  the  y  axis;  after  applying  this  roll,  OMEGA  is  the  angle  the 
x'  axis  pitches  up  from  the  x  axis;  positive  angular  deflections  are  defined  using  the 
Right  Hand  Rule.  These  angles  are  thoroughly  discussed  in  Section  2  of  the  text. 
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This  completes  the  type  E  component  input  data.  Now,  the  input  file  will  either  terminate  or 
the  next  entry  will  be  another  TYPE  variable,  followed  by  its  associated  component  information. 

Example  Input  format 

The  data  shown  below  is  a  listing  of  the  input  file  for  the  generic  submarine  described  in 
Section  10  of  the  text.  The  input  data  origin  is  implicitly  defined  as  the  nose  of  the  submarine 
by  XDIR,  the  setting  of  the  hull  shift  coordinates  to  zero,  and  the  first  line  of  the  hull  profile 
data.  Note  that  commenting  text  in  the  data  file,  other  than  the  COMMENT  input,  is  not  read 
since  ESAM  will  have  finished  reading  the  record  before  getting  to  it. 

8  #  of  components  making  up  the  submarine, 

n  Input  data  x-coordinate  increases  going  from  bow  to  stern, 

c  Output  data  x-axis  origin  is  opposite  centroid  of  ellipsoid  1 . 


0.0 

Output  data  z-axis  origin  is  coincident  with  input  origin 

453. 

Hull  length  used  for  nondimensionalization. 

hp 

HULL 

Component  1  is 

a  hull  component  of  primary  ’order*. 

/ 

shift 

coordinates  are  set  to  zero. 

23 

#  sections  ■> 

24  stations. 

0. 

0. 

0. 

0. 

1. 

6.7 

8. 

0.02 

2. 

0. 

11.2 

.04 

4. 

12. 

16.1 

.08 

8.5 

17.4 

21.5 

.17 

17. 

22.7 

27.4 

.35 

25. 

26.2 

31.2 

.51 

36. 

30.7 

33.7 

.72 

• 

73. 

38.1 

30.7 

1.5 

07. 

30.6 

41.4 

1.25 

102. 

40. 

42.1 

.0 

108. 

40.4 

44.2 

0. 

168. 

42.6 

40.0 

-2.5 

267. 

40. 

40.6 

-2.7 

200. 

38.4 

40.4 

-2.8 

200. 

37.7 

46.0 

-2. 

310. 

36.4 

42.0 

-.5 

313. 

36. 

40.6 

.2 

317. 

36.6 

37.4 

1.7 

363. 

20.0 

32.0 

.75 

388. 

21.2 

24.0 

.6 

410. 

14.0 

18.2 

0. 

440. 

6.5 

6.8 

0. 

463. 

0. 

0. 

0. 

Is 

SAIL 

COMPONENT  2 

168. 

.0..-27.5 

- 

00. 

60., -5. 4  50 . ,26.1  3., 26.1  .3 

Is 

STBD  SAILPLN 

3 

120. 

,7.73,-40.1 

0. 

20.5,-2.02  15.5,15.  6.0,15.  .15 

Is 

PORT  SAILPLN 

4 

120. 

,-7.73,-40. 

1  180. 

20.6,-2.02  15.5,15.  6.0,15.  .15 
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Is  TOP  RUDDER  5 


440. .0. .-3.4 

Is  BOT  RUDDER 

-90. 

6 

17.65,4.02 

14.35,18. 

o 

00 

.15 

440. ,0. .3.4 

Is  STBD  STRNPLN 

90. 

7 

17.65,4.02 

14.35,18. 

o 

M- 

00 

.15 

443. ,3. .0. 

Is  PORT  STRNPLN 

0. 

8 

29.  ,4.1 

19. ,22.3 

8.  ,22.3 

.15 

443. .-3. ,0. 

180. 

29.  .4.1 

19. .22.3 

8. ,22.3 

.16 

B.2  Output  Format 

The  output  file  (ESAM.OUT)  containing  the  output  from  a  run  using  the  above  input  is  shown 
on  the  next  page.  It  uses  a  full  132  character  line. 

Interpreting  this  output  file  is,  for  the  most  put,  straightforward.  Except  for  the  listing 
of  XOFF  and  ZOFF  in  line  5,  all  output  is  specified  relative  to  the  standard  body  fixed  system 
of  axes.  The  section  giving  the  ‘EQUIVALENT  ELLIPSOID  REPLACEMENT  GEOMETRY’  may  need 
some  explanation.  Ignoring  the  last  line  and  the  lines  beginning  with  ‘Optmzd  comp.  . 
this  section  is  simply  listing  the  characteristics  of  the  replacement  ellipsoids  for  the  specified 
components.  These  characteristics  determine  the  added  masses  of  a  component,  neglecting 
interference  effects.  The  last  line  of  the  section  gives  the  errors  associated  with  fitting  the 
asymmetyric  hull  with  a  single  ellipsoid,  as  defined  in  Section  8  of  the  text  by  equations  85 
and,  for  EXbar,  equation  77;  this  line  is  not  printed  for  a  hull  component  which  is  designated 
axisymmetric,  nor  for  one  of  secondary  order. 

If  a  component  is  of  secondary  order,  then  interference  effects  are  determined  by  the 
velocity  field  around  the  replacement  ellipsoid  of  its  associated  primary  component.  However, 
if  the  primary  component  is  of  type  H,  this  replacement  ellipsoid  geometry  is  first  optimized 
to  reflect  more  accurately  the  local  hull  component  geometry  opposite  the  x  location  of  the 
secondary  ellipsoid’s  centroid.  It  is  this  primary  ellipsoid  optimized  geometry  that  the  lines 
beginning  ‘Optmzed  comp  . .  .’  are  listing;  only  changes  to  the  original  replacement  ellipsoid 
geometry  are  noted  in  this  row.  This  optimized  geometry  is  only  used  for  calculating  the 
interference  velocities  for  the  secondary  component  of  the  previous  line;  it  has  no  effect  on  the 
added  masses  of  any  other  component,  including  the  primary  component  itself.  If  the  primary 
component  is  of  type  E,  no  optimization  takes  place. 
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